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INTRODUCTION 


In  the  last  decade,  much  effort  has  been  devoted  to  the 


obtaining  of  accurate,  non- relativistic,  Hartree-Fock  (HF)  represen- 

(1-5) 

tations  of  atomic  wave  functions.  For  closed  shell  systems, 

these  consist  of  a single,  energy  optimized.  Slater  determinant  made 
up  of  symmetry  adapted  spin-orbitals.  For  most  open  shell  cases, 
more  than  one  determinant  is  necessary  in  order  to  satisfy  the  sym- 
metry requirements  of  the  L,  S multiplets. 

( S V ) 

The  HF  method  is  useful  for  a large  variety  of  purposes, 

such  as  to  obtain  a qualitative  and  quantitative  description  of  optical 

and  X-ray  spectra,  or  to  evaluate  certain  one-electron  properties 

{ electronic  density  distributions,  diamagnetic  susceptibilities,  atomic 

scattering  factors,  etc.  ).  In  addition,  the  HF  solutions  serve  as 

excellent  starting  points  for  improved  calculations. 

The  accurate  computation  of  quantities  like  transition  proba- 

(8  9 ) 

bilities,  ’ inelastic  scattering  differential  cross  sections,  bind- 


ing energies  of  molecules, 
(11, 12) 


(10) 


electron  affinities  of  neutral  sys- 
tems, '*  hyperfine  structure  constants,  etc.,  and  the  proper 

understanding  of  phenomena  related  to  these,  however,  require  the 
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use  of  wave  functions  beyond  the  HF  approximation.  In  most  of 
these  cases,  a good  calculation  provides  a test,  not  only  for  the 
experimental  methods,  but  also  for  the  theory  itself. 

The  two  standard  variational  methods  for  improving  the  HF 

solutions  are:  (i)  the  interparticle  method,  in  which  the  expansion  of 

the  wave  function  Y includes  the  interparticle  coordinates  r-- 

ij 

explicitly,  ^ ^ and  (ii)  the  configuration  interaction  (Cl)  method, 

in  which  Y is  expanded  in  terms  of  properly  symmetrized  products 
of  one-electron  functions.  ^ ^ Both  have  been  proved  useful  in 

atomic  and  molecular  calculations,  the  first  one,  because  of  its  poten- 
tial accuracy,  the  second  one  because  it  is  easy  to  handle,  to  inter- 
pret, and  to  extend  for  larger  systems. 

The  aspects  involved  in  the  determination  of  N-electron  Y 's 
where  N equals  oris  greater  than  four,  are  of  an  essentially  different 
nature  in  the  two  methods.  In  the  interparticle  formalism,  the  prob- 
lems arising  in  the  evaluation  of  the  integrals  and  in  the  obtaining  of 
appropriate  computer  codes  to  handle  the  complicated  indexing  in- 
volved in  the  calculation  of  the  Hamiltonian  matrix,  are  not  entirely 
, , (17,18) 

solved  yet.  On  the  contrary,  there  are  no  more  real  numer- 

ical or  programming  difficulties  in  non- relativi stic  Cl  calculations  , 

(19,20)  , , 

and,  so,  at  the  present  time  Cl  is  the  most  promising  alter- 
native. In  this  work,  we  obtain  and  analyze  an  accurate  Cl  wave 
function  for  the  ground  state  of  the  carbon  atom. 

Chapters  I through  IV  are  of  a general  nature;  they  concern 
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the  theory  behind  the  procedure  used  in  this  calculation,  which  is 
applicable  to  all  first  row  atoms  as  well:  we  describe  in  detail  the 
optimization  of  the  parameters  of  the  Slater  type  orbitals  (STO's), 
the  use  of  first  order  density  matrices  of  selected  Cl  wave  functions 
to  obtain  the  one -electron  basis,  the  introduction  of  the  whole  degen- 
erate space  of  each  configuration  into  the  Cl  expansion,  and  the 

means  to  partition  it  to  improve  the  convergence  of  that  expansion. 

3 

The  calculation  of  the  P state  of  the  carbon  atom  is  de- 
scribed in  Chapter  V.  In  Chapter  VI  the  results  are  discussed  and 
compared  with  other  works  in  the  literature. 


CHAPTER  I 


THE  CONFIGURATION  INTERACTION  EXPANSION 


1 . 1 The  Method  of  Configuration  Interaction 

Let  be  a complete  and  orthonormal  set  of  one -elec- 

* 

tron  functions  of  the  space -spin  coordinate  x=  ( "r,  ? )'.  For  an  N- 

electron  system,  for  every  selection  of  N indices  k,  < k-><  . . .<k 

1 c n 

one  can  define  a normalized  Slater  determinant  : 


D. 


x^)  = (N!)-2  det  ’^k^^2^-- 


^k  (^n) 

n 


It  can  be  easily  proved  ^ ^ that  every  normalizable  antisym- 

metric wave  function  can  be  expressed  as  a linear  combination  of 
Slater  determinants  built  up  from  a complete  basic  set  of  one -elec- 
tron functions.  In  the  method  of  configuration  interaction  (Cl),  the 
exact  wave  function  is  hence  written  as 


Y (x^,  x^. 


■*N>  " 


s . 
1 


Di  (x^,X2, 


c. 

N'  1 


The  coefficients  c^  are  determined  variationally,  as  those  which 
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5 


minimize  the  total  energy  E : 

E = < T I H It  > = z ij  Hij  cfcj/  Zij  cfcj  . 

where  H-  = (D^[h  ! Dj  ) is  the  matrix  element  of  the  total 

Hamiltonian  H between  the  basic  Slater  determinants  D.  and 

1 

Dj  . The  condition  for  the  energy  to  be  an  extremum,  6E=  0, 
leads  to  the  well-known  matrix  eigenvalue  equation 

( IHl  - E n)  <c=  0 , 

which  has  non-trivial  solutions  only  for  values  of  E which  are 
roots  of  the  secular  equation 

det  ( HI  - E 11)  = 0 . 

For  a complete  set  of  one-electron  functions,  this  equation 
is  of  infinite  dimension,  and  cannot  be  solved  unless  it  is  truncated. 
What  is  done  in  practice  is  to  choose  a good,  small  one-electron 
set  {$^(1);  i=  1,Z,...,M^N}.  The  finite  Cl  expansion  consisting 
of  all  possible  Slater  determinants  formed  with  > is  called  a 

full  Cl  e^:pansion.  Both  the  eigenfunctions  and  the  eigenvalues  of 
the  full  HI  matrix  are  invariant  under  non- singula,r  linear  transfor- 
mations (n.  s . 1.  t . ' s)  of  the  set  because,  in  this  case  (sufficient 
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condition),  such  transformations  are  equivalent  to  n.  s.l.t.'s  of  the 
set  of  N-electron  determinants.  It  is  possible  to  find  a n.s.l.  t. 
of  the  original  set,  which  optimizes  the  convergence  of  the  full  Cl: 
this  is  discussed  in  Chapter  II.  A practical  definition  of  optimum 
convergence  is  given  in  Section  3.4. 


1.2.  Projection  Operators 

For  an  atom,  the  non- relativistic  Hamiltonian  H may  be 
written  as 


H = 


S h.  + 2 ..  g.. 

^ 1 i<j  ®ij 


where 


hj=  -1  7; 


IJ 


_Z^ 

r,- 


the  mass  polarization  terms  have  been  neglected,  and  atomic  units 
(a.u.)  are  used  throughout. 

This  Hamiltonian  has  certain  symmetry  properties  which  may 
best  be  taken  advantage  of  by  using  them  to  transform  the  set  of  N- 
electron  functions  {D^  } into  a more  convenient  set  } of  sym- 

metry adapted  functions.  Such  a set  may  be  obtained  by  means  of 

(16,21,22)  . 

projection  operators;  it  is  much  smaller  than  i } 
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because  it  spans  only  that  part  of  the  Hilbert  space  with  a certain 
symmetry  type. 

The  Hamiltonian  H defined  above,  commutes  with  the  oper- 

2 2 

ators  L , S , L^,  (squares  and  z-components  of  the  total  orbital 
and  spin  angular  momenta  ),  and  with  IP,  the  parity  operator.  Let 
us  denote  by  A any  one  of  these  operators.  Then,  for  non-degen- 
erate states  Y , 


A Y 


meaning  that  the  total  wave  function  Y is  an  eigenfunction  of  A 
with  eigenvalue  The  corresponding  projection  operator 

is  defined  by  the  property 


$ 


or  0 


where  ^ is  a normalizable  arbitrary  function,  and  satisfies 


An  explicit  expression  for  is  given  by  Ldwdin  : 


= n(  A-\^)/()^.x^) 
l^k 


(1.1  ) 
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since  A commutes  with  H , so  does  o!^  . In  addition. 

k 


satisfies  the  relations: 


°k  = 


0^=  0 for  k?5l  . 


1.3.  Projected  Configurations 


The  one-electron  functions  are  taken  to  be  symmetry  adapted 

spin-orbitals  the  form 

£ s 


=Rp(r)Y^  (0,cp)5  (cr) 

Irngm^  it'  ' Img  ' > t / ^ g/ 


where  the  ® ^ )'s  are  normalized  spherical  harmonics, 


(23) 


the  § {Cfg)'s  are  the  spin  functions  a or  8 , and  the  R.^(r)'s 

^ C 


are  linear  combinations  of  STO's,  S.ij: 


J(£) 


^i£  = ? Sjg  a.^. 


Sj£  = Nj£  r(^j  + e-2jf-r 

= [(2Z^j)  ^/(2£+2nj  + 2)  ! 


th 


e Ri£'s  being  orthonormal  among  themselves  for  each  £ value. 
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Let  us  now  define  a configuration  as  an  ordered  set  of  N 
indices  (i,  £•)  associated  with  the  spin-orbitals  i|».o  . To  each 

ilingnig 

pair  of  (i,  £ ) indices  there  correspond  many  (in  fact,  48+  2)  different 

spin-orbitals.  Hence,  one  can  form  several  Slater  determinants 

D^a  belonging  to  each  configuration  K and,  since  the  spin-orbitals 

are  symmetry  adapted,  the  determinants  are  eigenfunctions  to  L , 

z 

S^,  and  parity.  In  addition,  it  is  possible  to  form  certain  linear 
combinations  of  determinants  belonging  to  the  same  configura- 

tion  K,  such  that  the  are  eigenfunctions  to  L^  and  : we 

shall  call  them  projected  configurations.  When  only  one  linearly 

I 

independent  can  be  obtained,  the  configuration  K is  said  to  be 

non-degenerate  in  the  corresponding  L,  S state.  Otherwise,  a 

superscript  (p)  is  introduced  to  label  the  various  projected  functions, 

K is  then  called  g^^-degenerate. 

The  projected  configurations  are  symmetry  adapted  N-electron 

functions;  besides,  it  is  convenient  to  orthonormalize  them.  The 

orthonormal  set  of  eigenfunctions  to  L^,  S^,  L , S , and 

z z 

IP,  can  then  be  used  as  a basis  in  which  to  expand  the  total  wave 
function  Y of  the  same  symmetry  type: 


Y = 2^  a„ 

Kp  K Kp 


(1.2) 


The  are  called  terms  of  the  Cl  expansion.  They  are  obtained 


by  applying  the  projection  operator 
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o 


L^,S^ 

L,S 


(1.3) 


upon  the  determinants  > and  then  orthogonalizing  successively 

the  resulting  functions.  They  can  be  expressed  as 


$(P)  - 
K 


O 


L , S 
L,S 


[s 

a=l 


D 1 = 

Ka  Ka  ^ 


K 

S 

a=l 


D c^P) 
Ka^Ka 


(1.4) 


where  nj^  is  the  number  of  determinants  which  span  configuration  K 


In  the  basis 


the  matrix  elements  of  the 

KL 


Hamiltonian  are: 


(21) 


H 


(pq)  - 

KL  “ 


P 

2 

a=l 


(1.5  ) 


where  we  have  used  the  "turn-over"  rule  and  Eq.  (1.4).  The  turn- 
over rule  is  assumed  to  be  applied  in  the  most  efficient  way,  i.e. 


aP 


The  integrals  I 
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are  evaluated  as  described  in  Reference  23.  The  calculation  of 
the  and  c^J^^  coefficients  is  discussed  in  Reference  20. 

l\u  J\u 


1 . 4.  Classification  of  the  Configurations 


The  introduction  of  symmetry  adapted  spin-orbitals  permits 
a classification  of  configurations  into  classes  and  subclasses.  A 
class  is  defined  by  the  particular  partition  of  N among  the  spin- 
orbitals  of  different  £ values.  Within  each  class,  a subclass  is 
defined  by  the  particular  set  of  partitions  of  the  quantum  number  i 
among  the  spin- orbitals  with  the  same  t value.  Thus,  the  sub- 
classes (s^)^s^(p)^d  , s^s^Sy(p)^d  , (s^,)^Sj^p^p^d  , and  s^s^SyP^p^-d  , 
form  the  class  (3s,2p,  Id).  Under  linear  transformations  of  the  R^^ 
functions,  a given  configuration  is  transformed  into  other  configura- 
tions of  the  same  class.  A full  Cl  expansion  over  a given  class 
is  invariant  under  n.s.l.t. 's  of  the  provided  the  whole  degen- 

erate space  of  each  configuration  is  considered.  This  is,  however, 
a sufficient  but  not  necessary  condition,  since  if  only  a few  couplings 
- but  the  same  for  all  configurations  of  a class  - of  spin  and  orbital 
angular  momenta  are  introduced  in  the  Cl  expansion,  the  latter  is 
also  invariant  under  a general  n.s.l.t.  of  the  Finally,  some 

special  n.s.l.t.'s  can  be  found,  under  which  certain  truncated  Cl  expan- 
sions remain  invariant;  these  are  used  in  the  selection  of  the  one- 
electron  basis  ( see  Chapter  II),  with  the  purpose  of  improving  the 


rate  of  convergence  of  Y . 


CHAPTER  n 


THE  BASIC  ONE-ELECTRON  FUNCTIONS 

Any  complete  set  of  one-electron  functions  can  be  used  as  a 
basis  for  the  Cl  expansion,  and,  as  long  as  the  latter  extends  over 
the  complete  set,  the  exact  Y shall  be  obtained.  In  practice,  how- 
ever, it  is  important  to  obtain  a set  such  that  only  a few  terms 
provide  a relatively  accurate  ¥ . In  this  calculation,  a great  deal 
of  care  has  been  given  to  the  selection  of  the  basic  spin-orbitals. 

The  HF  approximation  for  closed  shell  systems  consists  of 
only  one  Slater  determinant.  When  it  is  used  as  the  zeroth-order 
wave  function,  all  the  other  determinants  in  the  Cl  expansion  can 
be  classified  as  one-,  two-,...,  n-excited  determinants,  according 
to  how  many  spin-orbitals  differ  from  those  in  the  HF  solution. 

For  non-closed  shell  systems,  where  the  HF  wave  function  may  be 
a linear  combination  of  determinants,  each  new  determinant  is  clas- 
sified as  above  with  respect  to  the  HF  determinant  it  differs  less 
from.  Let  us  point  out  here  the  difference  between  n-excited  deter- 
minants and  n-excited  configurations  ( or  n-excitations  ),  the  latter  being 
defined  as  configurations  differing  from  the  HF  one  by  n sets  of 
( i,  f ) indices . 
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It  is  also  customary  to  define  the  correlation  energy  E 

— — corr  ' 

as  the  difference  between  the  exact  non-relativistic  energy  E and 

nr 

the  HE  energy 

Ili? 

E = E - E„^ 
corr  nr  HE 

According  to  the  many-electron  theory,  the  HE  wave 

function  is  a good  starting  point  for  better  approximations , because 
it  has  the  formal  role  of  making  the  rest  of  the  exact  Y and  E 
come  out  simple:  the  Cl  expansion  consists  then  mainly,  of  HE  and 
two-excited  determinants,  while  the  correlation  energy  is  approxi- 
mated as  the  sum  of  the  correlation  energies  of  all  the  possible 
pairs  of  electrons.  In  Section  2.  4,  we  concentrate  on  obtaining  a 
fast  convergence  of  the  expansion  in  two-excited  determinants. 

Some  one -excited  ones  are  quite  important,  but  they  are  very  few, 
and  hence  do  not  constitute  a practical  problem.  The . convergence 
of  the  remainder  of  the  Cl  series,  which  consists  essentially  of 
three-  and  four-excited  determinants,  is  directly  related  to  the 
above,  as  will  be  discussed  in  Section  5.6. 

2.1.  Initial  Set  of  One-Electron  Eunctions 

We  start  then  with  HE  spin-orbitals  from  a nominal  HE  wave 
r • (25)  . , . 

function,  1.  e.  one  which  is  intermediate  in  accuracy  between  the 

minimal  STO  set  and  the  saturated  HE  wave  functions,  but  which  gives 
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an  energy  closer  to  that  obtained  from  the  latter  one.  For  the  first 
row  atoms,  the  nominal  HF  set  consists  of  a small  (4s,  3p)  STO 
basis.  Thus,  four  s-type  and  three  p-type  orbitals  can  be  formed 
from  it:  the  HF  orbitals  are  known;  the  others  are  found  by  succes- 
sive Schmidt  orthonormalization  of  linear  combinations  of  STO's.  A 
Cl  calculation  including  all  the  single  and  double  excitations  which  can 
be  formed  from  this  one-electron  basis  yields  a preliminary  Cl  wave 
function. 

The  second  step  is  to  increase  the  STO  basis.  The  experience 
of  many  worker s points  to  the  fact  that,  when  an  HF  solution 

is  used  as  the  zeroth-order  wave  function,  the  pair  correlations  are 
localized  largely  within  a shell  and  between  electrons  in  the  same  im- 
mediate region  of  space.  It  seems  hence  reasonable  to  look  for  new 
STO's  giving  rise  to  orbitals  which  are  optimum  for  expanding  the 
double  excitations  of  a pair  of  HF  spin-orbitals.  However,  this  re- 

V 

quires  the  obtaining  of  a large  number  of  STO's:  for  example,  one  of 

each  £- symmetry  type  for  every  pair;  furthermore,  many  of  them  have 

large  approximate  linear  dependencies  among  themselves,  and  this 

introduces  appreciable  errors  in  the  calculation  of  the  integrals.  What 

is  done,  then,  is  to  optimize  several  STO's  per  shell,  first  in  the  L 

shell  because  it  is  known  to  be  the  most  difficult  one,  then  in  the 

K shell  and  Intershell  alternatively;  within  each  of  these,  the  STO'S 

are  considered  in  order  of  expected  importance.  In  order  to  select 

a new  STO  one  chooses,  out  of  various  sets  of  parameters  (n.;Z.fl), 

J 
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the  one  which  yields  the  best  Cl  energy.  For  optimization  purposes, 
the  Cl  expansion  for  each  shell  includes  single  and  double  excitations 
alone,  and,  in  the  case  of  degenerate  configurations,  only  the  energy 
relevant  terms,  as  deduced  from  preliminary  calculations. 

Within  each  shell  and  symmetry  type,  "saturation"  with  STO's 
is  achieved  when  no  new  STO  can  be  found,  which  gives  any  substan- 
tial improvement  in  the  energy  (this,  of  course,  depends  on  the  par- 
ticular accuracy  desired). 

Let  us  assume,  then,  that  a satisfactory  STO  basis  has  been 
obtained.  The  corresponding  one-electron  basis  can  still  be  improved 
by  introducing  the  notion  of  natural  spin-orbitals. 

2.2.  The  Reduced  First  Order  Density  Matrix 

The  reduced  first  order  density  matrix  (1-ma.trix)  enters  our 
problem  when  we  try  to  improve  the  basic  spin-orbitals.  It  is  defined 
as(l^) 

Y(x^'|xj)=  N f ^(xj^x^.  . .Xj^)  dx^dx^.  . . dxj^  . 

For  a Cl  wave  Lmction  'f  = 2.  D.  c.  , where  the  determinants  D.  are 

111  1 

made  up  of  spin-orbitals 

Y(x'Jxp=  2 [2(-l)^ij  c^c  ] ^J(x')l(x  ) 
kl  ij 
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where  the  second  sum  goes  over  all  the  pairs  (ij)  of  determinants 
and  differing,  at  most,  in  the  spin-orbitals  k and  1 . This 

expression  may  be  written  as: 


The  matrix  ^ of  the  y,  , 's  can  be  easily  calculated.  Since  the 

kl 

generating  function  Y is  an  eigenfunction  of  IP  , L and  S , the 

z z 

operator  Y(l'|l)  commutes  with  IP{1),  L (1)  and  S (1);  this  follows 

z z 

immediately  from  the  properties  of  matrix  elements  between  Slater 
(23) 

determinants'.  The  matrix  ^ is  then  block-diagonal  in  the  follow- 

ing manner: 


oC  spin,  even  parity 


-1 

o 

-1 

0 

o 

1 

z 

> oC  spin,  odd  parity 


-3 


^ i 

: o 


O ! /3  spin 
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The  wave  function  Y is  assumed  to  be  an  eigenfunction  of 

2 2 

L and  S . If,  in  addition,  Y belongs  to  a state  with  M = 0,  the 

Lj 

block  associated  with  m£  is  identical  to  the  one  with  -m£,  for  a 

given  and  parity,  because  in  this  case,  Y is  invariant  under 

space  conjugation.  Analogously,  if  Y belongs  to  a state  with  M = 0, 

o 

the  block  with  CC  spin  is  equal  to  the  one  with  P spin.  In  con- 
nection with  the  selection  of  the  one-electron  basic  set  (see  Sections 

2.  3 and  5.  3),  these  properties  can  alv/ays  be  taken  advantage  of  for 

(23) 

the  following  reason:  it  can  be  proved'  ' that  the  Hamiltonian  matrix, 

2 2 

in  a basis  of  eigenfunctions  to  L and  S , is  independent  of  and 

Mg.  This  implies  that,  if  the  elements  of  degenerate  spaces,  for 

different  M values,  are  related  among  themselves  by 
Lj 


M^)  = [(L+ Mi^)(L±M^+  1)  ] ^ M^^il)  (2.1) 


and  an  analogous  relationship  holds  for  the  different  M values,  the 

S 

coefficients  a in  Eq.  (1.2)  are  the  same,  independently  of  M,-  and 
Kp 

M . In  Eq.  (2.1  ),  the  symbol  L,  stands  for  either  the  step-up  or  the 

o ^ 

step-down  orbital  angular  momentum  operator.  Therefore,  under  the 
conditions  above,  the  study  of  the  convergence  of  the  Cl  expansion  for 
an  L,  S state  can  be  performed  without  reference  to  the  M,  and  M_ 
value  s. 

When  Y belongs  to  a state  with  L=  0,  the  corresponding  1 -ma 

2 

trix  commutes  with  L (1):  then,  the  s and  d orbitals  do  not  mix. 
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neither  do  the  p's  and  f's,  and  all  the  blocks  of  a given  ^-symmetry 
type  are  equal. 

The  1 -matrix  can  be  diagonalized  by  a linear  transformation. 
Its  eigenvalues  n^^  are  called  occupation  numbers.  Its  eigenfunctions 
are  called  natural  spin-orbitals  (NSO's)  of  the  wave  function  from 
which  Y is  constructed.  The  1 -matrix  becomes 

Y(l'|l)=S  n.  xf(r)Xi(l) 

i 

2.  3.  Natural  Radial  Orbitals 

In  this  section,  we  consider  1 -matrices  derived  from  Y 's 
with  0 and  Mg=  0 (for  N odd,  it  is  not  possible  to  have  Mg=  0, 

but  the  considerations  below  can  still  be  applied  to  different  blocks). 

As  mentioned  in  Section  2.  2,  the  blocks  of  ^ of  the  same 
t- symmetry  type  are  equal  among  themselves,  if  Y belongs  to  an 
L=  0 state.  When  this  is  not  the  case,  ^ has  different  blocks  for 
different  values  of  ni^,  and,  hence , various  sets  of  radial  parts  of 
NSO's  spanning  the  space  defined  by  the  radial  parts  of  the  original 
spin- orbitals . The  question  is;  should  one  use  different  radial  func- 

tions to  represent  the  different  m£' s of  the  basic  spin-orbitals,  or 
should  one  choose,  perhaps  by  a maximum  density  criterion,  one  of 
the  sets  of  NSO's,  and  use  it  for  all  values  of  ni£? 

For  clarity,  let  us  select  a specific  ex3.m-ple:  we  consider  the 
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case  of  the  p-type  NSO's  of  a P wave  function  for  a two-electron 

atom.  There  are  two  sets  of  these:  the  R^j^'s,  corresponding  to 

m^il,  and  the  corresponding  to  m^^O,  where  the  subscript  i 

tV» 

labels  the  i ^ natural  spin-orbital  for  both  sets.  Let  us  form  the 
orbitals : 


pJ-RilYii  and  P?-RilY]_o. 


and  from  them,  a determinant  (p?p^/),  where  the  notation  is  such 
that  (i)  orbitals  to  the  left  of  the  bar  are  associated  with  OC  spin, 
those  to  the  right,  with  P spin,  and  (ii)  the  orbitals  are  positioned 
in  order  of  increasing  i , and,  within  a value  of  the  subscript  i , 
according  to  increasing  m^  values.  Let  us  now  obtain  the  '^P  pro- 
jection of  the  above  determinant: 


(pfPi‘/)=  (p?pj/)  - (p?pj/) 


where  it  has  been  necessary  to  introduce  the  following  functions: 


Pi'^ilYio  ^~^ilYii 


in  which  the  original  m^-  0 radial  part  is  mixed  with  the  m^=  1 
spherical  harmonic,  and  vice  versa.  Therefore,  unnecessary  redun- 
dancies  arise  in  the  description  of  eigenfunctions  to  L and  S"  when 
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the  true  NSO's  are  employed.  In  this  work,  we  have  chosen,  for 
each  f-value,  the  radial  functions  associated  with  the  set  with  largest 

leading  occupation  number,  and  used  it  for  all  rn^  values.  Following 

(28 ) 

Bunge,  ' ' these  one-electron  functions  shall  be  called  natural  radial 

orbitals  (NRO's). 

2.4.  Natural  Spin- Orbitals  for  Systems  with  more  than  Two  Electrons 

Lbwdin  has  suggested  that  the  NSO  expansions  of  many- 

electron  Y 's  might  have  optimum  convergence  properties  with  re- 
spect to  (i)  truncations  of  the  one-electron  basic  set,  and  (ii)  trun- 
cations of  the  Cl  expansion.  Both  of  these  properties  have  been 
proved  to  be  rigorously  true  for  two-electron  systems.  Concern- 

ing the  first  point,  (i),  our  experience  and  that  of  other  authors 
shows  that  a non-optimized  basis  may  be  drastically  truncated  by 
using  NSO's,  but  that  this  is  not  so  when  one  employs  the  spin-orbit- 
als formed  from  an  optimum  set  of  STO's,  as  is  the  case  in  the 
present  calculation.  As  for  (ii),  it  holds  practically  true  for  the 
ground  states  of  the  Be  and  Ne  atoms, ' which  have  symmetry, 
and  we  shall  try  to  find  a way  of  extending  it  to  other  symmetries 
as  well. 

Let  us  consider  a variational  wave  function  ^ , made  up  of 
the  HF  configuration  and  all  double  excitations  of  a pair  of  HF  spin- 
orbitals  , i'j  : 
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N 


N 


0(1,  2,.  . . ,N)=det  [n  l(i)}  + det{[(n  i)/t,  (k)l(l).]S  t (k)t  (l)c  ] 

i=l  ^ i=l  ^ ^ ^ n^N  P q Pq 


p>N 

q>N 


= 


(2.2) 


where  {t.;  i-  1 , . . . , N } spans  the  HF  space,  [ t ; p=N+  1,  . . . , M 1 
^ P 

defines  an  orthogonal  complement  to  the  HF  space,  and  the  Cp^'s 
are  variational  parameters.  The  corresponding  1 -matrix,  Y(0), 
has  the  same  eigenfunctions  and  eigenvalues  as  the  1 -matrix  derived 
from  ^(1,  2) 


cp(l,  2)=  det  [ t 1 } + det  { S t 't'  c } 

k 1 p q pq 


P>N 

q>N 


= det[tj^^^}+  X 


(k.l) 


where  the  Cp^'s  are  those  from  Eq.(2.2).  In  addition,  Y(0)  has 

the  eigenfunctions  ['1'^  ; i = 1,  . . . , N;  i^k;  i^l]  , each  with  eigenvalue- 

equal  to  one.  When  9(1,2)  is  expressed  in  terms  of  its  NSO's  (let 

us  call  them  (k,  1)  pair  NSO's),  is  transformed  into 

nat 

which,  according  to  the  two-electron  theory,  has  optimum  convergence 
properties  with  respect  to 

Next,  we  shall  consider  a new  variational  wave  function  ^ 

$ — ^ ^ ^(k>l)  _j_ 

HF 
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where,  for  definiteness,  k < 1<  m<  n^  N,  and  is  defined  in  a 

way  entirely  analogous  as  . Its  NSO's  shall  be  called  (k,  l)(m,  n) 

pair  NSO's.  Under  certain  conditions,  these  may  be  subdivided  in  two 
sets:  one  of  them,  similar  to  the  set  of  (k,  1)  pair  NSO's;  the  other, 
similar  to  the  set  of  (m,  n)  pair  NSO's.  This  implies  that  the 
(k,l)(m,  n)  pair  NSO's  are  nearly  optimum  for  expanding  both  pairs 
of  HF  spin-orbitals,  while  still  constituting  an  orthogonal  set.  It  has 
been  shown  that  such  a property  is  very  closely  fulfilled  for  the 
ground  state  of  Be, for  which  the  KL  shell  NSO's  are  obtained 
from  the  wave  function 

(1s)2(2s)2  ^ det{[2  ^S(x.  x.)k. . ](2s ) ^ } 

i s j ^ ^ 

+ det  [(ls)^[s  ^S(x.x.)L.  } 

1 J ij 

1 

One  can  expect  an  analogous  behaviour  for  the  joint  NSO's  of  any 
number  of  weakly  interacting  electron  pairs. 

Let  us  advance  the  hypothesis  that  the  NSO  expansions  have 
optimum  convergence  properties,  in  general,  for  any  group  of 
electrons  localized  in  the  same  region  of  space.  K this  is  true,  the 
joint  NSO's  of  weakly  interacting  groups  may  also  be  subdivided  in 
two  sets,  resembling  those  of  the  individual  groups.  For  atoms  in 
the  first  row  of  the  periodic  table,  there  are  two  distinct  groups, 
corresponding  to  the  K and  L shells. 
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The  alternatives  for  the  choice  of  the  basic  one -electron 
functions  are  as  follows; 

a)  K shell  NSO's;  optimum  convergence  in  the  K shell;  poor  con- 
vergence in  the  L shell. 

b)  L shell  NSO's;  poor  convergence  in  the  K shell;  optimum  con- 
vergence in  the  L shell. 

c)  KL  shell  NSO's;  good  convergence  for  both  shells. 

Since  the  triple  and  quadruple  excitations  appear  largely  as  products 
of  double  and  single  (or  double  and  double)  excitations,  a poor  con- 
vergence in  either  shell  means  a poor  convergence  of  these  as  well. 
As  for  the  Intershell  configurations,  any  basis  which  is  good  for  them 
appears  to  be  extremely  bad  everywhere  else,  and  they  shall  hence 
not  enter  our  discussion  at  this  point. 

In  the  next  section,  we  shall  describe  the  particular  choice  of 
NSO's  used  in  the  present  calculation. 

2.5.  Final  Set  of  One-Electron  Functions 

We  have  learned  by  experience  that  most  of  the  difficulties 
lie  in  the  L shell.  Thus,  in  a first  attempt,  we  concentrate  our 
^ efforts  in  obtaining  the  best  possible  set  of  one-electron  functions  in 
which  to  expand  the  L shell,  even  if  this  has  to  be  done  at  some 
expense  of  the  K shell  and  of  several  triple  and  quadruple  excitations. 
A KL  shell  NSO  basis  should  be  tried  next. 


Let  us  consider  all  the  orbitals  formed  from  the  HF  STO's 
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and  those  which  were  optimized  for  the  L shell.  A Cl  expansion  is 
then  set  up,  which  includes  almost  all  of  the  L shell  double  excita- 
tions, and  the  corresponding  L shell  NSO's  are  obtained.  (In  order 
to  preserve  the  HF  spin-orbitals  as  NSO's,  a few  configurations  are 
omitted:  this  is  illustrated  in  Chapter  V for  the  case  of  the  state 
of  the  carbon  atom.)  It  should  be  pointed  out  here  that,  except  for 
states,  one  does  not  deal  with  NSO's  but  with  NRO's,  as  explained 
in  Section  2.  3. 

The  L shell  NRO's  are  optimum  to  expand  the  L shell;  one 
can  only  hope  that  the  convergence  of  the  K shell  and  Intershell  ex- 
pansions in  that  basis  is  not  too  slow.  In  order  to  improve  the 
effects  of  those  configurations  containing  the  remaining  orbitals 
(i.e.  those  formed  from  STO's  optimized  in  the  K shell  and  Inter- 
shell), we  obtain  the  tail  NRO' s : these  are  the  NRO's  of  a Cl  wave 
function  consisting  solely  of  the  HF  configuration  plus  the  K shell 
and  Intershell  configurations  containing  the  "remaining  orbitals" 
described  above.  The  tail  NRO's  are  automatically  orthogonal  to 
the  L shell  ones,  since  the  two  sets  are  made  up  of  orbitals  which 


are  orthogonal  among  themselves. 


CHAPTER  m 


CRITERIA  FOR  THE  SELECTION  OF  CONFIGURATIONS 

Three  main  criteria  are  commonly  used  in  most  Cl  calcula- 
tions for  selecting  the  configurations,  and  the  choice  of  a criterion 
depends,  to  a large  extent,  on  the  particular  use  one  wishes  to 
make  of  the  wave  function,  and  also  on  the  personal  preferences 
of  each  author.  In  this  work  we  have  taken  into  account  all  three 
of  them,  at  various  stages  of  our  calculation.  Let  us  briefly  de- 
scribe them  here. 

3.1.  Second  Order  Perturbation  Theory 

In  the  Rayleigh-Schrddinger  (RS)  perturbation  theory,  the 
second  order  correction  to  the  energy  is  given  by  the  formula: 

E^^^=  2 H^./(E  -H..) 

oi  ' o ll' 

1 

where  H .=  (<!’  fH  'J'.  ),  ^ is  the  wave  function  for  the  unperturbed 
oi  o 1 o ^ 

Hamiltonian  H 

o 

H $ = E 
o o o o 
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and  defines  a complete  set  of  symmetry  adapted  functions 

spanning  the  orthogonal  complement  to  *^’o.  If  we  choose 

O fir 

and  [c2.  ,the  contribution  of  each  term  to  the  second  order 

1 ' K K 

energy  is  given  by 


This  expression  has  the  advantage  that  one  can  calculate  it  a priori , 
and,  hence,  have  an  idea  of  the  importance  of  the  HF  interacting 
terms  of  the  Cl  expansion. 

For  states  of  two-electron  systems,  the  RS  perturbation 

expansion  converges  fast,  and  in  addition,  the  e(3^'s  are 

Kp 

reasonably  good  approximations  to  the  energy  contributions  of  the 
individual  terms,  as  found  by  straight  application  of  the  variational 
method.  Furthermore,  their  goodness  improves  with  increasing 
nuclear  charge  Z. 

However,  for  outer  electron  pair  excitations,  the  RS  expan- 

(33)  (2) 

Sion  converges  slov/ly:'  ' E'  ' is  a poor  energy  correction,  and 

the  E^^^'s  are  often  not  even  of  the  correct  order  of  magnitude. 

Kp 

The  second  order  perturbation  theory  criterion  is  hence 
quite  rough,  and  it  can  only  be  used  in  the  very  first  stage  of  the 


calculation. 
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3.  2.  Eigenvector  Components 

Only  within  the  expansion  of  a given  pair  of  electrons,  does 
the  magnitude  of  the  eigenvector  component  provide  an  energy  cri- 
terion for  selecting  configurations:  as  pointed  out  in  Reference  20, 
these  magnitudes  are  larger  for  configurations  which  correlate 
outer  electrons,  than  for  those  correlating  electrons  closer  to  the 
nucleus,  when  the  energy  contributions  are  the  same. 

However,  if  one  is  interested  not  only  in  a good  energy, 
but  also  in  a good  wave  function,  terms  with  relatively  large  eigen- 
vector components  but  small  energy  contributions  should  be  pre- 
served, since  they  might  be  important  to  account  for  properties 
other  than  the  energy. 

3.3.  Partial  Energies 

Let  IHI  be  the  matrix  representation  of  the  Hamiltonian  in 
an  orthonormal  set  ; i=0, 1 , . . . } . The  eigenvalue  equation  is 

IHI  (c  = E (c  (3.1) 

If  a unit  row  vector  u=.100.  . . 0.  . . is  applied  from  the  left 
to  both  terms  of  Eq.  (3.1),  one  gets 

E=  S.  < 1h  I > c./c  = S.  E. 

1 O 1 1 O XI 
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The  's  defined  in  this  way  are  called  partial  energies  of  ^ 

ilC 

/ Q C \ 

with  respect  to  . ’ They  represent  the  actual  energy  contribu- 

tion of  ^ if  (sufficient  condition),  when  ‘I’  is  excluded  from  [$.  } 

K K.  X 

the  new  eigenvector  components  c!  satisfy 


c. 

1 


c. 

1 


for  all  i^  k 


(3.2) 


The  proof  is  very  simple.  Let  us  consider  two  secular  equations 
differing  only  in  that,  in  one  of  them,  ^ is  missing.  Then,  we 

iC 

have 


E(-$  ) =2  < <!>  |h  I $.  > c.'  / 


i 7^  k 


o '1  1 o 


= 2 <«l>  1h  1^.  )c./c  + <*  1h  |<l>,  > c,  /c 

o 'i  lo  o''k  k o 


E-E(-0^)=  = Ej^ 


The  fulfillment  of  Eq.(3.2  ) is  directly  related  to  the  stability  of  Y 

with  respect  to  the  exclusion  of  the  term  $ ; this  depends  on  the 

k 

nature  of  alone.  The  stability  of  y is  particularly  good  when 

. or,  even  better,  when  a normalized  truncation  of 

o HE  o TR’ 

a Cl  wave  function  which  includes  the  HE  configuration;  it  is  still 


further  improved  when  the  's  are  constructed  from  NRO's. 
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3.  4.  Optimum  Convergence 

In  our  criterion  for  optimum  convergence  of  a fixed  size  Cl 
expansion,  we  shall  take  into  consideration  both  the  eigenvector  com- 
ponent and  the  partial  energy  of  each  term. 

It  is  possible  to  find  a value  of  the  modulus  of  the  eigenvector 
component,  EIG,  below  which  the  terms  crowd  together,  most  of 
them  contributing  very  little  to  the  energy;  a few  of  these,  however, 
do  have  sizable  partial  energies,  and  one  can  set  a tolerance,  TOL, 
such  that  the  terms  for  which  the  absolute  value  of  the  partial  energy 
exceeds  TOL  are  separated  from  the  rest. 

In  our  234-terms  Cl  expansion  for  the  C atom,  we  shall 
preserve  all  the  terms  with  eigenvector  component  larger  than 
EIG=  0.0035,  and,  out  of  the  remaining  ones,  those  having  a partial 
energy  (in  absolute  value)  larger  than  TOL=0.  000012  a.u. 

The  convergence  is  arbitrarily  said  to  be  optimum,  if  no 
2 34-terms  Cl  wave  function  can  be  found,  which  gives  the  same  (or 
a better)  total  energy  for  smaller  EIG  and  TOL  values. 


CHAPTER  IV 


THE  BASIC  N-ELECTRON  FUNCTIONS 

As  mentioned  in  Section  1.3,  the  basic  N-electron  functions 
are  orthonormalized  projected  configurations,  or  terms,  . We 

shall  describe  here  how  to  form  them.  Except  for  the  determination 
of  degeneracies,  which  can  easily  be  done  by  hand,  the  obtaining  of 
the  terms  is  implemented  in  a computer  program  (see  Section  5.1). 

4.  1 . Degeneracy  of  a Configuration 

The  degeneracy  of  a configuration  K for  the  state  under  con- 
sideration determines  the  number  of  LS  terms  belonging  to  K . 

It  is  evaluated  in  the  following  manner;  first,  all  equivalent  spin- 
orbitals  are  combined  according  to  the  vector  coupling  model  and 
the  Pauli  principle, to  form  all  the  possible  multiplets.  Then, 
all  the  multiplets  arising  from  all  the  groups  of  equivalent  electrons 
are  in  turn  coupled  together,  in  any  order;  this  gives  a set  of  mml- 
tiplets,  some  of  them  appearing  several  times;  the  number  of  times 
a given  multiplet  can  be  formed  in  this  way,  is  the  degeneracy  g 
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of  the  configuration  K,  for  that  multiplet. 

Let  us  illustrate  this  procedure  with  a specific  example:  the 
degeneracy  of  the  configuration  (s^  s^s^  (p^  of  the  ^ state  of 
the  carbon  atom.  The  two  equivalent  spin- orbitals  ^^.n  only 

give  a S.  Hence,  the  must  be  formed  by  coupling  the  re- 
maining four  spin-orbitals:  (p^)^  gives  a , a and  a 

gives  a % and  a . The  two  possible  couplings  that  yield  a 
are: 


^S(s^^  ^^(3283)  ^(p^)^}  and  ^(8283)  V(p^)^ 


Thus,  the  degeneracy  is  equal  to  two. 


4.  2.  Projection  of  a Configuration 


In  the  case  of  the  orbital  and  spin  angular  momenta,  the 

projection  operator  of  Eq.(l.l)  can  be  written  in  a more 

A. 

(22) 

practical  form,  in  terms  of  the  well-know'n  step-up  and  step- 

down  operators  M,  and  M , defined  by  M , = M +i  M and 

+ - ■^  + X y 


M_=  M^-  i My.: 


T k -k 

Ayr  max 

(2k+l)(k+m)  ! 2 

’ /T,  \ t ^=0 

(k-m)  ! 


(-l)^M 


(k-m+V)^(k-m+v) 


(4.1) 


V ! (2k+v+l)  ! 


where  M stands  for  either  L or  S , k(k+l)  is  the  eigenvalue  of 
2 

M we  wish  to  project,  and  m is  the  eigenvalue  of  , 


which  is 
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assumed  to  be  known.  The  only  possible  k-values  range,  at  most, 

from  k=  m,  m+1,...  up  to  a certain  , which  may  be  evaluated 

by  the  addition  theorem  of  angular  momenta  and  which  depends  on 

the  particular  configuration. 

2 2 

The  operator  , defined  by  Eqs.(1.3)  and  (4.1),  can 

be  easily  applied  to  any  Slater  determinant  D : 

Kp 


. 2 D d<P> 

K S ^Kp  °Ka  ^Ka 


(4.2) 


In  Eq.  ( 4.  2 ),  n^^  is  the  number  of  determinants  which  span  the 

space  {K;  L,  S,  Mg  3;  not  all  the  D^^'s  with  the  correct  and 


Ka 


L 


i(p)  I, 


values,  belong  to  this  space;  also,  some  of  the  d'^'  's  may 

^ KcC 

equal  to  zero. 


be 


Let  us  assume  that  the  configuration  K is  degenerate,  and 
let  us  project  another  determinant  . If  the  corresponding 

projection  is  linearly  independent  of  the  two  functions 

may  be  Schmidt  orthogonalized;  if  not,  one  follows  the  same  proce- 
dure with  another  determinant,  and  so  on,  until  g independent 

^i(p)ts  have  been  found.  The  orthonormalized  functions  are  the 
x\ 

's  of  Eq.  ( 1. 4)  . 

K 

The  obtaining  of  the  independent  is  not  quite  as  "trial 

and  error"  as  might  be  deduced  from  the  above  description.  In 
fact,  one  soon  discovers  many  empirical  rules,  and,  with  a little 
practice,  one  guesses  immediately  which  determinants  to  project. 
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4.  3.  Partition  of  Degenerate  Spaces 


The  degenerate  configurations  would  not  create  any  special 
problems  were  it  not  for  the  fact  that  their  degeneracies,  in  the 
case  of  atoms  larger  than  Be  , are  often  very  large  (for  C , the 
largest  degeneracy  among  subclasses  giving  rise  to  double  excita- 
tions, is  27);  the  number  of  terms  to  be  considered  becomes  ex- 
tremely large,  and  a great  many  of  them  contribute  enough  to  the 
energy  so  that  they  may  not  be  ignored;  as  a result,  the  conver- 
gence of  the  Cl  expansion  is  exceedingly  slow. 

For  the  purpose  of  illustrating  this  problem  and  how  it  is 


carbon  atom,  which  is  6-degenerate  in  the  multiplet.  We  find 


solved  here,  let  us  consider  the  subclass  (s  )^(p  )^(d  of  the 

X ^y  z 


six  Slater  determinants  with  M =1  and  M =1,  whose  projections 
are  linearly  independent,  and  we  obtain  the  corresponding 

X\ 

When  these  are  introduced  into  the  secular  equation,  we  observe, 
in  general,  that  the  contribution  of  the  configuration  (s^^(p^  )^(d 
is  very  much  spread  over  all  of  the  terms.  Since  there  are 
6!  =720  different  ways  of  partitioning  the  6-degenerate  space,  just 
by  changing  the  sequence  in  which  the  projected  determinants  are 
successively  orthogonalized,  and  furthermore,  the  set  of  six  deter- 
minants is  not  unique,  one  speculates  that  some  choice  might  be 
such  that  only  one  or  two  of  the  terms  are  important. 

A given  ordered  and  complete  set  of  determinants  to  be 
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projected  for  a given  subclass  determines  a particular  partition  of 
the  degenerate  space  of  that  subclass.  We  shall  try  to  find  an  opti- 
mum partition. 

According  to  the  many-electron  theory,  one  expects  that 

terms  with  vanishing  Hamiltonian  matrix  elements  with  the  HF  con- 
figuration give  negligible  or  small  contributions  both  to  the  energy 
and  to  the  wave  function.  It  seems  reasonable  to  try  a partition 
in  which  the  first  five  determinants  to  be  projected  are  three-  and 
four -excited  determinants,  and  the  sixth  one  is  a two-excited  one. 

In  the  case  under  consideration,  for  example,  we  could  have: 


(sjpJdjd^/s^pJ) 

(Sjp|p|dj/Sid2) 

(s,pjd»dj/sip;) 
(Sjp[d|d“/a^pJ) 
(s  p^d^d^/s  p^) 

' ri 1 1 rr 

(s  p^p^d^/s  d^) 

1 r 1 1 1 1 


three  excited 

three 

four 

three 

three 

two 


It  should  be  noted  at  this  point,  that  it  is  not  always  possible  to 
partition  the  degenerate  space  in  this  way:  for  some  double  excita- 
tions, it  might  be  necessary  to  have  more  than  one  two-excited  de- 
terrhinant,  in  order  to  span  the  full  degenerate  space. 

Obviously,  the  partial  energy  contributions  for  the  first  five 
projections  (which  are  projections  of  linear  combinations  of  the  first 
five  determinants)  are  zero  since 
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<«=j^|h|^|p))=0  ; p=l,2,3,4,5. 


Now,  the  interaction  of  a configuration  with  the  HF  configuration  is 
a constant  given  by: 


g 


K 


^F 


where  aj^p  is  the  coefficient  of  the  HF  configuration,  and  the  a 's 

are  those  of  the  terms  of  configuration  K in  the  final  wave 

function.  This  implies  that  must  have  maximum  interaction 

K 

with  the  HF  configuration,  provided  that,  under  any  other  partitions, 
the  E^^'s  be  always  negative  or  zero,  which  is  the  case  for  most 


configurations.  Of  course,  the  partial  energies  with  respect  to 
(see  Section  3.3),  will  not  be  necessarily  zero  for  the  first  five 
terms:  in  the  example  shown  above,  we  have  only  succeeded  in 
finding  an  appropriate  partition  in  which  all  the  two-body  correla- 
tions (two -excited  determinants)  are  lumped  together  into  the  last 


TR 


element  of  the  degenerate  configuration. 


CHAPTER  V 


CALCULATION 

5.1.  The  Computer  Program 

Our  Cl  calculation  of  the  ground  state  of  the  carbon  atom  is 
performed  with  a computer  program  originally  written  by  Bunge, 
and  later  improved  with  the  addition  of  a projection  operator  section 
and  modified  for  use  with  highly  degenerate  configurations.  Let  us 
summarize  briefly  the  main  parts  of  the  program  in  its  present 
form. 

There  is  first  a very  simple  input,  which  specifies  the  type 
of  work  requested  (Cl  + configuration  search;  optimization  of  STO 
parameters;  Cl  + corre  spending  ^ and  NRO's;  ...,  etc.),  the  number 
of  electrons  N , the  maximiim  size  of  the  secular  equation  during 
search,  the  LS  state,  the  STO  parameters  and  coefficients,  the  deter- 
minants to  be  projected,  the  configurations  and  their  degeneracies, 
and  finally,  the  terms  to  be  inserted  in  the  Cl  expansion.  Unless 
previously  calculated  and  stored  on  magnetic  tape,  the  one-  and  two- 
electron  integrals  are  evaluated  in  double  precision  arithmetic.  The 
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next  step  is  to  generate  and/or  read  the  "Master  Tape",  which  con- 
tains all  the  relevant  information  about  the  subclasses,  their  b and 
c coefficients  and  corresponding  Slater  determinants,  and  all  the 
configurations  considered.  To  each  term  the  program  associates  a 
number,  which  is  all  one  needs  to  specify  in  the  input  when  such 
term  is  wanted  in  the  Cl  expansion.  Any  selection  of  terms  is  per- 
missible. The  Hamiltonian  matrix  is  computed  as  described  in  Sec- 
tion 1.3,  and  diagonalized  by  the  Givens  method.  Finally,  if 

desired,  the  1 -matrix  is  calculated  and  the  NRO's  are  obtained.  If 
instead,  one  wishes  to  search  for  the  most  important  terms  out  of 
a large  number  of  them,  the  program  can  discard  certain  terms, 
according  to  the  tolerance  in  the  eigenvector  components  and  the 
partial  energies  indicated  in  the  input,  pick  up  new  terms,  obtain 
the  corresponding  E and  Y , discard  some  terms  again,  . . . , and 
so  on  until  all  the  terms  have  been  tried. 

The  last  part  of  our  calculation  was  performed  on  a CDC 
3600/3400  electronic  computer,  with  41  K available  words.  The 
program,  as  adapted  to  this  machine,  permits  a maximum  size  of 
the  secular  equation  of  234x234,  a maximum  degeneracy  of  27,  and 
no  more  than  400  Slater  determinants  in  the  expansion  of  the  pro- 
jection of  a given  determinant.  In  addition,  for  each  configuration, 
the  number  of  terms  inserted  times  the  size  of  the  secular  equation, 
may  not  exceed  4000.  The  lowest  eigenvalue  and  its  corresponding 
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eigenvector  can  be  found  in  less  than  seven  minutes,  while  the  gener- 
ation of  the  Hamiltonian  matrix  takes  between  three  and  forty  minutes, 
depending  on  the  configurations  involved. 

5.  2.  The  STO  Basis 

We  use  as  a starting  point,  the  nominal  HF  wave  function 

(25) 

calculated  by  Bagus,  ' ' with  an  STO  basis  consisting  of  four  s- 

and  three  p-type  STO's,  denoted  as  follows: 


for  £=  0 , SjQ  with  j = 1,2,3,  and  4 , 

for  1 , Sj  with  j=  1,2,  and  3. 


In  addition  to  the  three  HF  orbitals  s^^,  s^,  Pj^ , four  more  orbitals 
can  be  formed  from  such  a basis,  i.e.  s , s , p , and  p . We 
arbitrarily  postulate  an  s^  consisting  only  of  S^q.  and  orthogonalize 
it  to  the  HF  orbitals:  s^  becomes  a linear  combination  of  all  four 
S^q's.  The  remaining  orbitals  are  similarly  obtained. 

Some  preliminary  L shell  NRO's  are  then  constructed  from 
the  (4s,  3p)  orbital  basis,  in  the  way  outlined  in  Section  2.  5,  thereby 
removing  the  arbitrariness  in  the  orthogonalizations  above  (the  L 
shell  NRO's  are  uniquely  defined,  since  they  are  derived  from  a 
wave  function  which  is  invariant  under  linear  transformations  of 
the  symmetry  spin-orbitals  other  than  the  HF  ones).  A full  L shell 
Cl  calculation  with  respect  to  single  and  double  excitations  is  carried 


out  with  these  "preliminary"  L shell  NRO's,  and  a selection. 


39 

L (4,  3,0,0),  is  made  of  the  most  important  terms;  the  letter  L in- 
dicates that  the  terms  are  L shell  excitations,  and  the  ciphers 
inside  the  parentheses  denote,  in  order,  the  number  of  s,  p,  d,  and 
f orbitals  involved. 

Let  us  now  energy  optimize  a d-type  orbital,  d^  , consisting 
of  a single  STO,  , using  an  expansion  which  includes  L(4,  3,  0,  0) 

plus  all  the  L shell  single  and  double  excitations  that  can  be  formed 
with  d^  . Next,  the  expansion  L(4,  3,1,0)  is  used  as  a starting  point 
to  optimize  , and  in  a similar  way,  the  process  is  applied  suc- 
cessively to  > and  . Any  further  attempts  to  ex- 

tend the  L shell  STO  basis,  give  negligible  improvements  in  the 
energy;  the  "saturation”  remainders  are  about  0.00010  a.  u.  for  s- 
type  STO's,  0.00010  a.u.,  0.00017  a.u.  , and  0.00010  a.u.  for  p- 
d- , and  f-type  STO's  respectively. 

It  is  interesting  to  note  that,  often,  more  than  one  absolute 
minimum  can  be  obtained  for  E,  as  a function  of  the  STO  parameters. 
Let  us  denote  the  latter  by  ( m ; Zj£  ),  and  let  us  consider  the  opti- 
mization of  S^.  We  find  three  sets  (1  ; 1.3),  ( 1 ; 2.  3),  and  (2  ; 2.4) 
located  in  slightly  different  regions  along  the  r axis,  yet  giving  the 
same  energy  minimum.  We  choose  (2  ; 2.4)  because  it  is  more 
orthogonal  to  the  other  S^^  functions,  and  therefore  one  can  expect 
a greater  flexibility  in  the  Cl  expansion. 

The  final  L shell  NRO's  are  now  evaluated  as  described  in  the 
next  section:  from  here  on,  the  symbols  s ^ , . . . , s p , . . . , p^,  d^  , and 
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<^2  denote  L shell  NRO's.  The  subscripts  are  used  to  classify  the 
NRO's  according  to  decreasing  occupation  numbers  of  the  1 -matrix 
block  from  which  they  are  selected.  The  f-type  NRO's  need  not 
be  obtained  because  very  few  terms  containing  f^  and  f^  appear  in 
the  Cl  expansion. 

Next,  let  us  consider  the  Intershell  expansion  (s,s_  and 

1.  w 

®1  ?!  excitations).  Two  parallel  Cl  calculations  are  performed, 

one  involving  Intershell  excitations  alone,  1(5,4,  2,2),  and  one  on 

both  L shell  and  Intershell  excitations,  LI(5,  4,  2,  2).  A comparison 

of  the  partial  energies  of  the  Intershell  terms  in  both  Cl  wave 

functions,  shows  no  appreciable  differences,  and  hence,  1(5,  4,  2,  2) 

is  used  as  a starting  point  for  the  optimizations  of  , and  S 

60  51 

(the  necessity  for  the  superscript  I,  meaning  Intershell,  will  be 
apparent  below). 

Independently,  we  consider  the  K shell.  A comparison  of 
the  partial  energies  of  the  K shell  terms  in  K(5,4,  2,2)  and 
KL(5,4,  2,  2),  justifies  our  using  K(5,4,  2,  2)  for  the  optimizations 
of  S^K),  s»),  sm.  S(K),  S(K),  sTO,  S(K).  Once  si^). 
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60 


and  have  been  found,  the  K shell  NRO's  are  obtained: 

in  this  entirely  new  orbital  basis,  K(6,  6,  2,  2)  converges  very  fast, 
w'hich  simplifies  enormously  the  search  for  the  remaining  S.^g^^'s. 
The  K shell  NRO's,  however,  are  no  good  for  expanding  the  L 
shell,  and  hence  they  will  not  be  used  any  further. 

When  ‘ ^61^^  used  in  an  Intershell 
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Cl,  the  results  are  even  better  than  those  obtained  with  and 

These  STO's  are  hence  abandoned  in  favor  of  their  K shell 
counterparts,  and  the  superscripts  are  dropped.  Also,  no  or 

are  considered,  assuming  that  their  K shell  counterparts  do 
an  equally  good  job. 

The  overall  "saturation”  remainders  for  the  Inter  shell  are 
probably  smaller  than  0.00020  a.u.  For  the  K shell,  they  can  be 
estimated  very  accurately,  from  the  spdf  energy  limits, 

and  the  computation  of  the  exclusion  effects  (see  Chapter  VI). 

The  simultaneous  optimization  of  two  STO's  is  found  to  be 
useful  only  in  the  case  of  ^2Z  ’ gs-in  over 

the  corresponding  successive  optimization  is  0.000050  a.u. 

The  STO  parameters  are  given  in  Table  1. 


5.  3.  The  NRO  Basis 


In  order  to  obtain  a set  of  L shell  NRO's  among  which  the 
three  HF  orbitals  appear  unchanged,  we  select  a Cl  expansion 
containing  the  HF  configuration  and  the  L shell  terms,  but  from 
which  certain  subclasses  and  configurations  have  been  banished. 
There  must  be  no  two  configurations  differing  in  only  one  orbital, 
where  the  latter  is  either  one  of  the  HF  orbitals.  For  example, 
the  configuration  (s  ) s s (p  ) interacts  with  the  HF  configuration 

X w X X 

2 2 2 

(s  ) (s  ) (p,  ) to  produce  a mixing  of  s_  with  s.,  i.e.  the  second 
12  1 2 X 

2 3 

s-type  NRO  is  a linear  combination  of  s^  and  s^ ; also,  (Sj^)  (pj^)  p^ 
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TABLE  1 
STO  BASIS 


j 

n£ 

Shell 

1 . 

1 s 

9.  153 

K (HF) 

2 

1 s 

5.  382 

K (HF) 

3 

3s 

3.  076 

L (HF) 

4 

2s 

1.428 

L (HF) 

5 

3 s 

2.  200 

L 

6 

2s 

10. 500 

K 

7 

3 s 

12. 000 

K 

1 

2p 

5.  152 

L (HF) 

2 

2p 

2.  177 

L (HF) 

3 

2p 

1.150 

L (HF) 

4 

3p 

2.  300 

L 

5 

2p 

9.  000 

K 

6 

3p 

11. 000 

K 

1 

3d 

1. 950 

1. 

2 

4d 

2.  000 

L 

3 

3d 

11.800 

K 

4 

4d 

13. 000 

K 

1 

4f 

2.  500 

L 

2 

5f 

4.  100 

L 

3 

4f 

16. 000 

K 

a.  The  usual  (n£)  notation  is  employed;  where 

b.  Localized  mainly  in  the  shell  indicated. 


n=  n.  + {+  1 , 
J 
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2 4 

interacts  with  (s^  (p^ ) to  mix  and  p^. 

In  addition,  as  pointed  out  in  Section  2.  2,  there  shall  be  a 

mixing  of  s-  and  d-type  orbitals,  and  of  p-  and  f-type  orbitals, 

unless  the  state  under  consideration  has  L=  0.  This  is  due  to 

2 2 

interactions  between  configurations  such  as  (s.)  s d (p,  ) and 

1 X y ^1 

(s,  )^(s  )^(p,  )^ ; however,  the  most  important  term  of  (s,)^s  d (p, 
1X1  lx  y 1 

can  be  written  as  a linear  combination  of  determinants  which  are 

2 2 2 

two-excited  with  respect  to  (s^)  (s^)  (p^)  : the  s^'^^y  could 

be  avoided,  were  it  not  for  the  presence  of  (s^  )^(d^)^(p^)^  and 
2 2 

(Sf)  ■ Although  some  of  the  configurations  of  the  type 

2 2 

(Sf)  Sx^y(Pj^)  extremely  important,  we  prefer  to  obtain  the 

NRO's  which  optimize  the  convergence  of  the  Cl  expansion  with 
respect  to  the  highly  degenerate  ( s )^(d^)^(p^  and  (s^  )^d^d^(p^  )^  . 

For  the  calculation  of  the  matrix  ^ , it  is  convenient,  as 

explained  in  Sections  2.  2 and  2.  3,  to  use  projected  configurations 
corresponding  to  0 and  0.  Of  course,  this  involves  a 

great  deal  of  additional  work,  inasmuch  as  one  has  to  select  a whole 
new  set  of  projected  determinants  for  each  subclass.  However,  % 
is  greatly  simplified,  since  only  the  block  with  OC  spin  needs  to  be 
calculated,  and  within  it,  for  each  parity,  only  those  blocks  with 
different  mg  values.  In  our  case,  since  we  have  used  only  s,  p, 
and  d orbitals  in  the  calculation  of  ^ , there  are  six  different 

blocks  instead  of  eighteen.  The  main  advantage  • of  this  manner 
of  proceeding,  is  that  it  makes  a given  choice  of  NRO's  optimum 
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for  a larger  number  of  spin-orbitals;  for  example,  for  the  d-type 
orbitals,  there  are  three  different  blocks  corresponding  to  |m.  1=  % 
and  Im^jsO,!,  and  2,  and  hence  three  (instead  of  ten)  different  sets 
of  radial  parts  of  NSO's  from  which  to  select  the  set  with  largest 
leading  occupation  number.  For  the  p-type  NRO's,  the  |m^|=l 
block  is  used;  for  the  d-type  ones,  the  [m^j^Z  block. 

The  tail  NRO's  of  the  K shell  are  obtained  by  using  a "tail" 
Cl  expansion,  consisting  solely  of  the  HF  configuration  and  all  the 
K shell  configurations  formed  from  the  following  basic  orbitals:  s^  , 

. s?  , Pi  , p^  , p^  , d^  , and  d^,  with  the  exception  of  the  HF  mixing 
and  the  s-d  mixing  configurations.  Of  course,  the  energy  of  such 
a Cl  wave  function  is  expected  to  be  quite  bad.  The  advantage  in 
obtaining  the  tail  NRO's  is,  that  the  convergence  of  the  final  Cl 
expansion  is  improved  with  respect  to  these  orbitals,  when  the  pro- 
portion in  which  the  tail  configurations  enter  the  final  Cl  wave  func- 
tion is  the  same  as  in  the  tail  Cl.  For  instance,  let  us  consider 


the  K shell  excitations  (d^)^,  (d^)^,  and  d^d^ : the  only  relevant 
terms  are  those  for  which  the  K shell  remains  coupled  to  symme- 
try, and  thus  the  analogy  with  two-electron  systems  is  valid,  i.  e. 
the  contribution  of  d^d^-in  terms  of  tail  NRO's  vanishes  identically, 
while  that  of  (d^)^  is  increased.  The  tail  NRO's  would  probably  be 
more  interesting  if  s^,  s^,and  p^  were  allowed  to  participate  in  the 
tail  Cl  expansion. 

The  STO  coefficients  of  the  final  basis  are  given  in  Table  2. 


ORBITAL  EXPANSION  COEFFICIENTS  IN  TERMS  OF  STO  BASIS 
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5.4.  The  Terms 


Let  us  consider  the  very  important  3-degenerate  configuration 
2 2 

(Sf)  complete  determinantal  space  of  this  configu- 


ration, for  M^=l  and  Mg=l,  is  made  up  of  14  determinants;  it 
involves  both  one-,  two-,  and  three-excited  determinants.  As  ex- 
plained in  Section  4.  3,  one  looks  for  sets  of  three  linearly  independ- 
ent projected  determinants,  containing  as  many  three -excited  ones 
as  possible:  in  this  case,  only  one.  We  chose  the  set: 


No  special  reason  will  be  given  for  this  choice;  it  can  only  be  based 
on  trial  and  error,  and  it  would  take  too  long  to  find  the  optimum 
set,  out  of  all  the  possible  ones  having  a three-excited  determinant 
as  their  first  component. 

The  advantages  of  the  procedure  outlined  in  Section  4.3  for 
partitioning  the  degenerate  space,  become  much  more  obvious  when 
one  considers  configurations  such  as  ^2^^!  ^^^2^3  ’ 

degenerate.  It  is  possible,  in  this  case,  to  find  25  linearly  independ 

ent  projections  of  three-  and  four-excited  determinants.  The  26^^ 
tv» 

and  27  components  of  the  set  are  projections  of  two-excited  deter- 
minants, and  they  are  found  to  correspond  to  the  couplings 


two 


three  excited 


one 
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S (s  j S2)  % (P2P2)  } ^ (Pj  (Sj^  S2)  (p^p^)  ^(Pj)^  respec- 

tively. 

For  most  configurations  involving  K shell  excitations 

2 2 

(xy)(s2)  (Pj)  , the  terms  have  been  chosen  so  that  only  one  element 

makes  any  contribution  to  the  energy:  it  invariably  represents  a 

coupling  ^S(xy)  -^(pj^)^,  reflecting  the  strong  tendency  of 

the  K shell  to  maintain  the  symmetry.  Some  core  polarization 

is  effected  by  configurations  of  the  type  s^(s2)^  (p^  d and 

2 2 

(®2^  (Pf ) Px^v  ’ pairs  (s  d ) and  (p  f ) cannot  be  coupled 

y X y X y 

to  a . 

For  L,  shell  configurations,  the  contributions  of  three-  and 
four -excited  determinants  associated  with  double  excitations  are  not 
riegligible.  For  those  subclasses,  when  we  use  the  same  partition 
of  the  degenerate  space  for  both  shells,  we  find  that,  although  the 
term  with  coupling  %(xy)^(p^)^  contributes  the  largest 

partial  energy,  all  the  other  terms  are  often  non-negligible  too.  ' 

It  is  necessary  then,  to  look  for  an  additional  optimum  partition 
within  the  remainder  of  the  degenerate  space:  this  is  done,  mainly, 
by  changing  the  order  in  which  the  three-  and  four-  excited  deter- 
minants are  projected;  occasionally,  different  sets  of  projected 
determinants  are  tried.  In  general,  all  one  achieves  is  to  increase 
the  partial  energy  contribution  of  the  largest  terms  at  the  expense 
of  the  smaller  ones.  For  configurations  which  are  not  very  impor- 
tant, this  is  often  enough  to  permit  us  to  discard  most  of  the  terms. 
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The  (p^f  excitations  are  all  non-degenerate,  and  so  the 
1 21  2 3_ 

coupling  S(Sj^)  ^(s^)  ■^(xy)  is  the  only  possible  one. 

For  single  excitations  such  as  (s^  S2S^(Pj  f , the  dominant 

term  does  not  belong  to  the  coupling  ^S(s  ^S(s  s ) f 

2 3 1. 

but,  rather,  to  ^(s^)^  % (s^s^)  (p^  )^  } . This  can  be  ex- 

plained in  the  following  way:  the  first  coupling  involves  only  one- 
excited  determinants,  while  the  second  one  involves  both  one-  and 
two-excited  ones,  the  latter  causing  greater  interaction  with  the 
HF  configuration. 

The  terms  of  the  subclasses  appearing  in  the  final  Cl 
expansion  are  given  in  Tables  3,  4,  and  5.  The  subclasses  are 
characterized  by  letter  subscripts;  the  determinants  are  given  for 
the  simplest  configuration  of  each  subclass,  rather  than  in  general, 
because  it  seems  easier  to  read  them  that  way;  they  are  written 
in  the  notation  introduced  in  Section  2.  3 , and  they  are  listed  in 
the  order  in  which  they  are  projected.  Since  the  projection  coef- 
ficients are  the  result  of  a series  of  projections  and  ortho- 

normalizations, it  is  not  always  possible  to  write  them  in  the  form 
of  quotients  of  simple  integers  or  radicals;  in  such  cases,  they 
are  given  in  decimal  form,  the  number  of  figures  depending  on  the 
accuracy;  when  all  the  coefficients  of  the  projection  have  the  same 
denominator,  the  latter  is  written  only  once,  for  the  first  deter- 
minant. 


The  first  ten  subclasses  of  Table  3 are  common  to  the  K 
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TABLE  3 


SLATER  DETERMINANTS  AND  COEFFICIENTS  OF 

Kcc 


L AND  K SHELL  CONFIGURATIONS 


2 2 2 

Subclass  Ll=  K1  : (s^  (s^)  (p)  , 1 determinant,  degeneracy  = 1, 


Determinants 


1 / 0 1/  \ 

1 (SjSjPjPj/SjSj) 


Ka 


2 2 

Subclass  L2=  K2  : (s^)  s^s^(p)  , 4 determinants,  degeneracy=  2. 


i Determinants 

1 / 0/  ^ 

1 (S1S2S3P1/S1P1) 

->  ! 0 1/  \ 

2 (S^S^PjPj/SjSj) 


Ka 

2 


b<2) 

Ka 

2/'v/2~ 

2 


Z . A 
s 

Determinants  b' 


Subclass  L3~  K3  : (s)  (p)  , 1 determinant,  degeneracy=  1. 

,(1) 


1 (Sjp}pfpJ/s^pl) 


Ka 


1 


2 3 

Subclass  L4=  K4  : (s)  (p^)  determinants,  degeneracy  = 3. 


i Determinants 

1 (SipJp}p2/Sjp}) 

2 (s^p|pjp2/s^p}) 

3 ( = iPiP}p^/^P?) 


b(l) 

Ka 

5/n/To 


b(2) 

Ka 


Ka 


-1 

-2 

2 


/ 
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Subclass  L5=K5:  (s)^  (p^)^ 


TABLE  3 continued 

12  determinants,  degeneracy  = 5. 


i 

Determinants 

'^KO. 

b(5) 

Ka 

1 

( s p^  p^p^  / s p*^) 

0 

10/n/33 

2/n/3 

2 

(SlP?P^P^/sjP°) 

-1/nT2 

0 

0 

3 

< = ip1p2P^  = ip1> 

2 

11 

1 

4 

(Sipjp‘p‘/Sipj) 

22 

2 

5 

, 0 10/  Ox 

(-JPJPJP^/^IP^) 

3 

Subclass  L6=  K6  : (s 

>^<Pw>^PxPy  • 

28  determinants. 

degi 

i 

Determinants 

1 

(SiP°P2P3/sjp}) 

-6/n/m 

2 

(s,pJp/p^/Sjp“) 

-10 

3 

(Sip1pJp3/s,P^) 

-10 

4 

(Sipjp^p^/Sipj) 

- 6 

5 

^®1^1^2^3^‘‘’l^l^ 

- 3 

6 

(SjpJpPp^/s^p}) 

- 3 

7 

<^p1piPe'  = iP3> 

- 5 

8 

(SjpJp[p“/Sjp“) 

12 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinant:  (s j^p^p^p^ / s j^p^) 


TABLE  3 continued 


Subclass  L7-  K7  : (s^)^  iPyf'  > 27  determinants,  degeneracy=  6. 


i 

Determinants 

b<2) 

Ka 

b(3) 

Ka 

1 

0.  781508148 

0. 361963863 

2 

(Sip\p}dj/sjd^) 

2.  526215823 

-0. 518629654 

3 

(s.pjdjdj/.jp;) 

1. 930314651 

4a 

1 

-2. 923333333 

1. 801963268 

- 5/s[45 

2 

-0. 776059887 

-0. 724407157 

3 

-1. 112537987 

-0. 190117275 

0 

4 

(=,p;d|dO/s^pj) 

2. 750684472 

-2. 946817765 

0 

5 

(s^pjd^d^/s^pj) 

3. 796013990 

5 

6 

(s^pjp}d^/sjd^) 

15 

Subclass  L8=K8:  (s)^ 

(p)^dd  , 58 

determinants,  degeneracy  = 10 

i 

Determinants 

X y 

b(5) 

Ka 

b(8) 

Ka 

b(io) 

Ka 

1 

0.  2086440 

-3. 6547708 

\FT5/4iz 

2 

(sip\d|d^/s  jp}) 

-1. 0432202 

3. 5540488 

^TGO 

3 

, 0 ,0  ,1  , 0> 
(SiPjdid^/SjPi) 

0. 

5.  2163146 

0 

4 

(s^pJpJd^/sjdj) 

-0. 4818428 

-2. 8909798 

- \[M 

5 

, ljO,0/  0, 

(SjP^d^d^/SiPi) 

3. 5408069 

-0. 1220935 

^/l20 

6 

(SipJpJdpSid^) 

-2.  1295515 

- ^^60 

7 

(sipjd^dl/sipj) 

-1. 1654978 

n/TE 

52 


TABLE  3 continued 
Subclass  L8=K8  (continued). 


i 

Determinants 

'ia 

C 

8 

, 1^0,1/  L 

(SiPidid^/^iPi) 

6. 1679395 

sT^/’sUz 

9 

(SiPidid^/SjPj) 

-n/240 

10 

(s^pJpJdJ/sjd^) 

n/i20 

Subclass  L9=K9:  (s)^ (p)^  (f )^,  42 

determinants,  degeneracy  = 6. 

i 

Determinant  s 

'ia 

b(6) 

Ktt 

1 

(SiPjf|£^/sjpJ) 

^/^/^/T5 

0. 

0 

2 

(s  p^f^f^/s  p*^) 

' 1^1  1 r Kr 

0. 

0 

3 

(Sip\ i/sip}) 

4TW 

0. 

0 

4 

(Sipfplf/Zsifj") 

-1. 4819172 

0 

5 

(SjP^fJ/SifO) 

2.  8166174 

n/T4/n/T5 

6 

(s^p“p;fj/s^fj) 

^105 

Subclass  L10=K10: 

, 14  determinants. 

degeneracy=  3 

i 

Determinants 

1 

(SiS^pjd^/sipO) 

-n/3/n/3 

-^/^5/^/3~ 

2 

(s^s^pJpJ/sjdJ) 

sTso 

3 

4Jo 

TABLE  3 continued 


Subclass  LI  1 : (s  ) 
' u' 

s s p p , 8 

w x^y^z ’ 

determinants,  degeneracy  = 

i 

Determinants 

41' 

41' 

1 

(S1S2S3PO/S1P}) 

1AT3 

1 

2 

(s^  S3pJp^/s^S2) 

3 

1 

3 

<=l=2=3P^  = lP2> 

2 

Subclass  L12  : (s 
w 

SxPyP^-d  , 32  determinants,  degeneracy 

i 

Determinants 

41' 

4t' 

1 

(SiS^pJd^/sipQ) 

2 

(s^pjp^dj/s^s^) 

-\/l690 

3 

<'i  = 2P2'*I'  = ip1> 

- n/640 

4 

(SiS^plp^SidJ) 

n/4840 

4a 

1 

n/5/6 

^B/^Uz 

2 

-n/To 

3 

V40 

-n/40 

4 

n/1000 

n/160 

5 

( s s p^p^  / s d*^) 
' 1 2^1^2  1 L 

n/1080 

n/120 

6 

(s^plpldj/s^s^) 

n/i60 

TABLE  3 continued 


2 

Subclass  L13  : (s^)  s^pdf  , 56  determinants,  degeneracy=  6. 


i 

Determinants 

Ka 

Ka 

1 

^idj) 

1. 059148182 

n/TT5’/6 

2 

(^iS^dffJ/ 

= lP?) 

0. 599144689 

n/14/3 

3 

(SlP}d[£[/ 

"l=2> 

3. 555121506 

^/2625/18 

4 

^®1  ^2^1'^!  ^ 

.0, 

s,£i> 

n/^/3 

^a 

1 

n/7/n/T8' 

n/ZT/n/Ts 

2 

^^^/3 

n/42 

3 

^/35 

4 

^/4^ 

5 

(SlSjpJdJ/: 

SlfJ) 

n/TO/2 

n/T^ 

6 

(s^  S2  p^^f  °/: 

SidJ) 

Subclass  L14  : PyP^  ’ ^ determinants,  degeneracy  = 

i Determinants 


1 (SjSjpJp'/SjSj) 


2 2 2 

Subclass  LI  5 : (s  ) (s  ) (d)  , 2 determinants,  degeneracy^  1. 


w ' X' 

Determinants 


1 


Ka 

5/^/^5 


TABLE  3 continued 


Subclass  L16:  , 


3 determinants,  degeneracy=  1. 


i Determinants 


1 (Sj=^£»£j/SjS2)  7/-^l 


Subclass  L17:  (s)^(p)^f 
i Determinants 


(s,PiP,f  !/s 


iPl^ 


, 9 determinants,  degeneracy^  1. 


4T5 


Subclass  L18 : ( s , 44  determinants,  degeneracy=  4. 


Determinants 


1 (SjpJp°f  J/sjP°) 


-4‘c! 


35/n/T^ 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinants:  (s^  pj  p{  f s^  p^),  ( s ^pjp^f  J / s^  p j ), 
(=iPMp“/=ifj). 


Subclas  s LI  9 : (s^)  determinants,  degeneracy^  2. 

i Determinants  bi3  ^ 

Ka 

1 (s  S2d^dJ  / Sj^d^)  n/tO/3 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 

T 0 2 

the  following  determinant:  (sj^  d^dj^dj^ /s^  S2)  . 


TABLE  3 


continued 


Subclass  L20  : (s)  p^p^(d)  , 64  determinants,  degeneracy  = 11, 


i 

Determinants 

b(5) 

°Ka 

b(io) 

°Ka 

1 

(SjpJp^dJ/sj^dJ) 

1. 5683137 

0. 

2 

(s^pJpOdJ/s^d}) 

0.4527332 

1. 5249857 

3 

(Sipld0d2/sip0) 

2. 5699830 

1. 0791262 

4 

(Sipld}dO/sip}) 

1. 0393190 

0. 5576189 

5 

^®1^2*^l'^l  ^^^1^1^ 

3.4985711 

0. 8626622 

6 

-0. 9149914 

7 

(s^pjd^dj/sjpl) 

-0. 1078328 

8 

(Sip^d^d^/sip  °) 

-0. 7624929 

9 

(SiP>0dJ/Sid|) 

-0. 5282705 

10 

(Si  plpld}/s^d2) 

3. 2349832 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinant:  (Sj^p^^dj-dj- /sj^p®)  . 


Subclass  L21  : (s)^p(d)^f  , 118  determinants,  degeneracy^  12. 


i Determinants 
1 (s^p°djd}/s^f°) 


4. 3204938 


The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 

the  following  determinants:  (s^p^d^dj/s^f  ^),  (s^  Pj^ d J d J / s ^f  ^), 

' Ic  /|2^l4:3/^  „1^  ^2^1^2/^  ^0  1 o -:> 


TABLE  4 

SLATER  DETERMINANTS  AND  b(P)  COEFFICIENTS  OF 

KCt 

INTERSHELL  CONFIGURATIONS 


Subclass  II:  s„s„(s„)^{p)^  , same  as  L2. 
' w X y 


w X'  y'  ' 

Subclass  12:  s„;s^s„s^(p)^ 

, 15 

determinants,  degeneracy^  6. 

i 

Determinants 

b(5) 

Ka 

b(M 

Kd 

1 

(S1S2S3S4/PJP}) 

l/^^3 

-1 

2 

(SiS2S4p}/s3pO) 

-3 

1 

3 

(S1S2S3PJ/S4P}) 

-3 

3 

4 

0 

-2 

5 

(s^S2p{p}/s3S^) 

3 

-1 

6 

(s^S3pJp}/s2S^) 

0 

2 

Subclass  13:  s s (p)"^  , 4 
w 

determinants,  degeneracy=  2. 

i 

Determinants 

b(2) 

Ka 

1 

(sjS2p}p\/p°p}) 

2/^/2 

2 

(sip}pfp}/s2p}) 

2 

•5 

Subclass  14:  s s (p  ) p 
w x'^y' 

, 27 

determinants,  degeneracy  = 9 

i 

Determinants 

b(8) 

Ka 

b(9) 

Ka 

1 

(s^S2p}p2/p?pJ) 

4/nT6~  -4/n/2 

2 

(siS2p}pi/pJp2) 

6 

-2 

3 

(Sip}p}p^/S2p}) 

0 

6 
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TABLE  4 continued 


Subclass  14  (continued). 


i 

Determinants 

b(8) 

Ka 

b(9) 

Ka 

4 

(s^S2pJp}/p}p°) 

9/^i6 

3/slz 

5 

(BjS^PiP^/PiPi) 

5 

1 

6 

(SiS^PiP^/PiPi) 

10 

1 

7 

(S2p}pfp2/sipj) 

0 

4 

8 

( s s p^p^ /p ^p 

\ J 2^L1 

6 

2 

9 

(Si  pfp/p^Vs^p/) 

0 

4 

Subclass  15:  s^s^(Py)^(P2)^  > determinants,  degeneracy  = 13. 


i 

Determinants 

b(^) 

Ka 

1 

(Si  S2p/p2Vp°P2^) 

-Z/\l6 

2 

, 0 1/00, 
(SjS^PlP^/PjPp 

4 

3 

/ 0 0 1/  0, 
<"lPl_P2P2'=2Pl> 

-2 

4 

( Pj'pj  Pj  / ®2  ) 

-2 

5 

(=1  SjPlV/p/pj) 

6 

6 

/ 0 0 1/  0, 
(S2P1P2P2/S1  Pi) 

2 

7 

/ 0 10/  0, 
(SiPiPiP2/=2P2> 

6 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition. 


the  following  determinants:  (s^  s^pj  p° /p}  P^)-  (Sj  s^p^P^/p}  pj  )> 
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TABLE  4 continued 


Subclass  I6:  s s (p 

w x'^y' 

PzPu  ’ 

determinants, 

degeneracy  = 

i 

Determinants 

b<2) 

°Ka 

b(2M 

KO 

b(27) 

°Ka 

1 

(SlS2p}p‘/pJp*) 

-2/n/T^ 

0 

0 

2 

(^iPjp^Pj/s^pJ) 

21 

0 

6/^/3 

3 

(SjS^pJp^/pJp') 

0 

6 

4 

(s^pJpjPj/SiPj) 

-12/3 

- 

8 

5 

<=1=2P2P3'p1p?> 

13 

7 

6 

^®1  ®2P2P3'^^i^l^ 

21 

7 

7 

(SiS2P^P^/p°p{) 

- 6 

- 

6 

8 

(^2Pip;p2/=iPi» 

6 

- 

2 

9 

'®2P{P2P3^"lPf> 

-15 

- 

1 

10 

*=iP1'’2'’^=2'’1’ 

6 

4 

11 

(S1PJP1P3/S2P2) 

6 

0 

12 

, 0 0/  1 Ov 

<=l"2PlP3^PlP2> 

0 

6 

13 

(S^pjpjpl/Sipl) 

6 

0 

14 

<"lPlP2P3^"2Pl> 

-15 

- 

13 

15 

(s  p^p^p^/s  p^) 
' D1^2^3  2^r 

- 3 

- 

5 

16 

( s p ^p^p ^ / s p 

^ r r2^3  2^r 

- 3 

1 

17 

( s p®p^p^ / s p 

' DL2^3  2^1 

6 

2 

18 

( s s p^ p^ /p ^p^) 

' 1 2^D2  3' 

-12 

- 

2 

19 

(s  p^p°p°/s  p°) 
' 2^r  2^3  r 1 

3 

- 

1 

20 

(SiP^PiP2/s2P3^ 

9 

3 

21 

(SipjpiP2/s2P3^ 

3 

- 

1 

TABLE  4 continued 


Subclass  I6  (continued). 


i 

Determinants 

b(2) 

Ka  °Ka 

Ka 

22 

(s  p^p^p^/s  pM 
' 1^1^1^3  2^2^ 

15/3 

1/nT3 

23 

(s  s p*^p^/p®p^) 
1 2^L3  2 

12 

8 

24 

(SiS2pJp^/pJp°) 

- 6 

- 4 

25  ■ 

^®2^1^2^3^®lPl^ 

-27 

- 9 

26 

(s  s p^p^/p^p^) 
' 1 2 ^1  ^1  ^2  ^3  ^ 

9 

3 

27 

/ 0 10/  Ov 

'"lPlPlP2'  = 2f3> 

6 

Subclass  17:  s s (p) 
w 

^(d)^,  100  determinants, 

degeneracy 

= 15. 

i 

Determinants 

b'"» 

Ka 

1 

(SiS^pJdJ/pJdJ) 

^/30/^/3 

2 

(Sip}p}d0/S2d}) 

-n/TO 

3 

(SjS^pidJ/pJdJ) 

-n/30 

4 

(SjS2pId2/p0dJ, 

-n/20 

5 

6 

(siS2P|<i[/p;dj) 

-n/40 

7 

(SipJpJdJ/s^dJ) 

^/30 

The  full  degenerate 

space  is  obtained  by  projecting,  in  addition 

the 

following  determinants:  (Sj^  s 2P^d^/p Jd ^ ), 

/ , 1 _ 0 j 0 / 

s^d^). 

TABLE  4 continued 


Subclass 

18  : 

(s  r s s p p , 
u w x^y^z 

same 

as 

Lll 

Subclass 

19  : 

same 

as 

L12 

Subclass 

110 

same 

as 

L13 
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TABLE  5 


SLATER  DETERMINANTS  AND  COEFFICIENTS  OF 

KL  SHELL  CONFIGURATIONS 


Subclass  KLl  : (p^)  (Pjj)  » determinants,  degeneracy  = 5. 
i Determinants 


1 (pjpjpjp'/p;?*)  -i/-^3 

2 (pjp?p;p»/p;p“)  3 


The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinants:  (pjp^^p°p^/pypj ),  (pjpjp°p^/pjp J ), 

, 0 1 T 0/  0 L 

(PlPlP2P2/PiPi)- 


Subclass  KL2  : (p^)  P^Py  ’ determinants,  degeneracy  = 9. 

(j) 

1 Determinants  bLi 

CL 

1 (pjpjp^pj/pipj)  2/^/6 

2 (pypjp^p^/pjpj)  6 

, / T 1 1 0/  T L - 

3 (P1P1P2P3/P1P1)  2 

4 (pjpjp^p^/pyp})  2 

C , T 0 1 0/  0 L 

5 (P1P1P2P3/P1P1)  2 

6 (pj-pjpjp^/pjp^)  -2 

7 (p[p?p;P°/p1p^)  6 

The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinants:  (pjpjp^p^/pjp j ),  (p jpjpjp^/pjp^)- 


TABLE  5 continued 


Subclass  KL3  : (pj^  (d)^  , 


27  determinants,  degeneracy=  6. 


i  Determinants 


1 (pjpjp}dj/pjdj)  - 5/^ 

2 (Pi^]°p}dJ/pJdJ)  15 


The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 
the  following  determinants:  (p^pjdjd  J /p^pj  ),  (p^pjdjd  j /pjpj  ), 


(p°p^d^d°/p°p^), 


(p^p^d^d^  /p°P^)- 

'^1^1  1 r^l^L 


Subclass  KL4 : s s s (p)^d  , 
w X y'^' 


53  determinants. 


Determinants 

°Ka 

jSjSjd^/pJp") 

- 25/sJlWo 

- Z5/\f3b 

iSjpJpJ/s^d^) 

25/N/3000 

jSjpfdJ/s^pl) 

^320/3 

degeneracy=  9. 


The  full  degenerate  space  is  obtained  by  projecting,  in  addition. 


the  following  determinants:  (s^  s^s^pj /pjdj ),  ( s^  s^p^dj  / s^pj  ), 

(SiS2p}d2/s3pJ),  (s^SjpJdJ/s^pJ),  (s^s^pjpj/s^dj). 


Subclass  KL5  : (p  f (p  p^  , 85  determinants,  degeneracy  = 21 . 

i Determinants  bl3^^ 

Ka 

1 (p|p2P2P3/PiPi) 

2 (pJp{p2P3/p}p2) 

3 (p}p}p2P3/p}p2' 


TABLE  5 continued 


Subclass  KL  5 (continued). 


11 


13 


17 


18 


Determinants 


4 (p[pjp“p“/p[p^) 

5 (p1p|p»P°/pJp‘) 

/p1p“) 


(P1P1P2P3 


7 . (p|p1p^P?/p?P°) 

8 (plpJplpl/pJpO) 

9 (pjpjp^pj/pjp^) 

10  (P?P2P2P3/p}pJ) 


(P°P°P^P^/P}P?) 


19  / 1 1 0 1 / 1 1, 

12  (PjP^P2P2/PiP3) 


(pJp1p2P3/p1p2) 


14  (plpOplpO/pOpl) 

15  (p}pjplp0/p}pl) 

16  (P?p1pIp^/p?P°) 


(pipjpjp^/p'p*) 

(p^p^p^p^/p^p*^) 

1 r 2^3  r 2 


14/nT3 


- 2 


12 


- 8 


-12 


- 4 


- 3 


10 


The  full  degenerate  space  is  obtained  by  projecting,  in  addition, 


the  following  determinants:  (pJp^P^Ps/pfpJ  )>  {pfp^P^Ps/pf  p}  )> 

(P1P?P2P3/P}P2)- 
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TABLE  5 continued 


Subclass  KL6:  s (p^)^  ^ determinants,  degeneracy=  26. 

i Determinants  b^3  ^ b^^^  b^^^ 

Ka  Ka  Ka 


1 

(SipJp^/pJpO) 

0. 779965319 

0. 832276131 

2 

(SjpJp°d}/p}pl) 

0. 573125876 

0. 350645971 

3 

(SjpJpld}/p}pO) 

0. 573125876 

0. 350645971 

4 

(sip}p}p^/p^d}) 

2. 579066443 

1. 099556503 

5 

(Sjp}p0dJ/pJp0) 

^(8) 

^Ka 

.(11) 

^Ka 

2. 870102209 

1 

-0. 800249258 

0. 493336795 

2 

-0. 044509342 

-0. 071435955 

3 

-2. 149303254 

0. 027005928 

4 

0. 303762687 

0. 085555993 

5 

-0. 650920044 

-0. 920828417 

6 

(SiPiP^di/piP^) 

-2. 255568812 

-0. 766619174 

7 

^®lPlP2^r  ^1^2^ 

-3. 332158754 

-1. 187570110 

8 

(s  p^p^d°/p°p^) 
' 1^1^2  r^i^2^ 

3. 928165571 

-0. 077579517 

9 

(SipJpOdO/pJpl) 

-0.  851584885 

10 

^®1  ^1^2*^!  '^^1^2^ 

-0. 166382392 

11 

(s^p}pldj/p0pl) 

3. 487230624 

The  full  degenerate 

space  is  obtained 

by  projecting, 

in  addition. 

the  following  determinants:  (s p^p^d^/pj p^),  (s jpJpjp^/p^d^ ), 
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TABLE  5 continued 


Subclass  KL6  (continued). 


(Si?! P2^1 /?! P2)>  (s^pjp2^2 /pjp^)* 


(SipJpX/p/pz*' 

/ 0 1 0 j 2 / 0 \ 

<PlP2P2'‘l/=lPP' 

/ 0 10/  OjOx 

(SiPi  pjp^/p^d^). 


/ 1 0 ,1  / 0 1, 

(s^p^p^dj/PiPi), 

(pJpJp'dJ/Sjp'), 


(^iPjpJdf/pJp?). 


2^2' 


/ 10  1/0,1, 

(SiPiPiP2/P2di), 

1-.1 


/ 0 10/iji, 

(SjPiP2P2/Pidi), 


(P^P^P°P^/s 
^^1  M ^2^2  1 1 ’’ 


I 0 1 lj2,  1, 

<PlPlP2'^l/“lP2>' 


and  L shells  (they  are  denoted  by  LI  Kl,...,  LIO  KIO);  for  LI 
through  L9,  the  projection  corresponding  to  the  largest  value  of 
the  superscript  (p)  , always  belongs  to  the  coupling  ^S(s  f'  ^S(xy) 
The  other  projections  correspond  to  less  obvious,  or 
linear  combinations  of,  internal  couplings.  Table  4 contains  the 
subclasses  of  the  Intershell  excitations:  the  only  reported  element 
of  13,  15,  and  17,  corresponds  to  ^SCs^s^)  %3(p^)^,  while 

for  12,  14,  and  I6,  the  two  elements  in  the  table  belong  to 
{%(sj  S2)  %(xy)  } ^(p^)2  and  (s^  S2)  % (xy)  (p^  )2  , in  that 
order.  In  Table  5,  the  subclasses  from  which  only  KL  triple 
and  quadruple  ( and  larger)  excitations  can  be  obtained,  are  reported; 
for  KLl , KL2,  KL3,  and  KL5,  the  couplings  are  of  the  same  type 
as  in  the  HF  configuration:  *S  (p^)?  *S(pj)^  ^(Pj)^  . (p^p  ) *S  (pj 
V(pj)2  , ‘S(d)f  Is(pj)?  %>(pj)2  , and  ‘s  (p/  ‘s  (PjP^)  %>  (Pj  ^ , 
re  spectively. 

Only  about  a half  of  all  the  subclasses  for  which  the  best' 
possible  partition  of  the  degenerate  space  was  obtained,  appear  in 
the  final  Cl  expansion,  and  are  reported  here. 

5.5.  The  Search  for  the  Terms 

As  mentioned  in  Section  5.  2 , the  three  large  divisions  of 
the  Cl  expansion,  K shell,  L shell,  and  Intershell  excitations, 
are  quite  independent  of  each  other,  inasmuch  as  the  partial  energies 
and  eigenvector  components  of  one  of  them  are  not  appreciably 
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affected  by  the  presence  of  the  other  two.  The  K shell,  L shell, 
and  Intershell  Cl  expansions  are,  hence,  investigated  first  inde- 
pendently; the  best  energies  obtained  for  each  of  them  are  indicated 
in  Table  6 , together  with  the  energies  of  their  most  relevant  trun- 
cations. The  last  column  of  Table  6 gives  the  error  in  the  trunca- 
tions, as  estimated  from  the  search  by  adding  all  the  partial  ener- 
gies of  the  terms  which  were  discarded  (this  estimate  is  usually 
good  to  within  10  per  cent). 

Let  us  consider  the  L shell.  The  search  has  to  be  done  by 
parts,  according  to  the  limitations  of  the  computer  program,  as 
indicated  in  Section  5.  1.  At  this  early  stage,  it  is  not  necessary 
to  set  the  tolerances  in  the  eigenvector  components  and  the  partial 
energies,  very  low;  we  are  interested  mainly,  in  obtaining  a rough 
approximation  to  the  value  of  the  partial  energy  contribution  of 
each  term;  we  use  EIG=  0.0035  and  TOL=  0.  000080  a.  u. 

The  K shell  Cl  search  consists  of  approximately  720  terms. 
The  convergence  is  good,  and  the  final  K shell  Cl  expansion  con- 
tains 49  terms.  This  is  much  more  than  the  necessary  number  of 
terms  when  a K shell  NRO  basis  is  used;  yet,  it  does  not  compare 
so  unfavorably  with  the  K shell  Cl  expansion  obtained  for  Be,  with 
a KL  shell  NSO  basis:^^^^  the  main  disadvantage  lies  in  its  effect 
on  the  convergence  of  the  KL  triple  and  quadruple  excitations. 

The  search  of  the  Intershell  configurations  yields  a large 
number  of  terms,  whose  contributions  are  small  but  not  negligible. 
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The  search  is  slowed  down  by  the  occurrence  of  large  degeneracies; 

this  limits  the  size  of  the  secular  equation,  as  explained  in  Section 

5.  1:  for  example,  when  searching  the  27-degenerate  subclass, 

2 

®l®2^Pp  PxPy  ’ maximum  number  of  terms  which  can  be  intro- 
duced at  one  time  is  4000/27  , i.  e.  148.  The  total  truncation 
error  in  this  case,  has  been  computed  by  adding  the  following  par- 
tial energies: 


(S1S2) 

d d 
X y 

-. 000300 

(S1S2) 

(f)^.  f f 

' ' ’ X y 

-. 000100 

(SiS^)  - 

others 

-. 000025 

(s^Pi)  — 

all 

-.000025 

When  the  K shell  and  Intershell  expansions  are  put  together, 
the  results  are  remarkably  additive.  This  does  not  happen,  however, 
when  the  K and  L shell  expansions  are  combined:  in  this  case, 
although  the  partial  energies  of  the  K shell  excitations  do  not  change 
much,  those  of  the  L,  shell  ones  are  substantially  diminished.  The 
total  energy  loss  by  non-additivity  is  about  -.0025  a.u.  , and  it  is 
partially  recovered  when  the  triple  and  quadruple  excitations  involving 
both  the  K and  L shells,  are  introduced. 

Finally,  all  three  expansions,  and  the  triple  and  quadruple 
KL  excitations  are  lumped  together  into  a final  search,  to  obtain 
the  total  wave  function.  The  corresponding  truncation  error  reported 
in  Table  6,  is  calculated  by  adding  the  following  partial  energies: 
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K and  L shell  truncations 


-.000943 


Intershell  truncations 


-.000645 


Triple  and  quadruple  KL  truncations  -.000315 


Triple  and  quadruple  KL  from 
subclasses  not  considered 


-.000500 


5.6.  Results 


Our  final  Cl  wave  function  consists  of  234  terms  formed 
from  a determinantal  space  of  dimension  4836.  Its  corresponding 
energy  eigenvalue  is  E = -37.  83378  a.u. 

The  eigenvector  component  of  the  HF  configuration  is 
0.972691373.  The  rest  of  the  wave  function  is  displayed  in  Tables 
7 through  18. 

The  first  seven  of  these  tables  contain  excitations  of  the  L 
shell  alone.  All  the  HF  interacting  (S2)  excitations  appearing  in 
Table  7 belong  to  the  coupling  ( s ^ (xy)  ^(p^^.  The  HF  non- 
interacting excitations  are  significant  only  for  the  pairs  (82)^  and 
S2P^  • In  these  cases,  and  also  in  Tables  13  and  18,  the  partial 

energies  are  computed  with  respect  to  ^ ; all  the  others  are 

TR 

given  with  respect  to  ^ . The  last  six  terms  of  Table  12  have 


zero  or  positive  partial  energies;  the  corresponding  values  with 

respect  to  # are:  -.000083,  -.000102,  -.000043,  -.000032, 

1 K, 


HF 


-.000030,  and  -.000020  a.u.,  respectively. 
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TABLE  7 


EXPANSION  OF  HF  INTERACTING  (s^)^ 

EXCITATIONS 

C onfiguration^ 

Degenerate^ 

Element 

Eigenvector 

Partial  Energy 

. 130723807 

-15401 

1=^2 

3 

038569453 

- 2371 

6 

032093068 

- 3779 

(S3) 

030227981 

- 2759 

P1P3 

3 

. 022263542 

- 1370 

(P^)" 

5 

018990747 

- 1074 

6 

-.009298258 

- 750 

(d^)2 

6 

-. 005675494 

- 313 

-. 004830124 

- 177 

(P3)^ 

5 

-. 004305142 

- 167 

P2P3 

8 

. 003742989 

19 

P2P4 

8 

. 003022259 

- 151 

(f  / 

6 

-. 002914846 

- 149 

10 

-. 001767170 

- 110 

P2P5 

8 

-. 001006853 

14 

(P4)^ 

5 

-. 000976456 

38 

(33)^ 

-. 000478943 

22 

a.  Only  the  excited  (i,  £)  indices  are  shown. 

b.  The  ordering  of  the  degenerate  elements  is  given  in  Table  3. 
c . In  units  of  10"^  a.u.(C). 
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TABLE  8 


EXPANSION 

OF  HF  NON 

-INTERACTING 

(s^)^  EXCITATIONS 

Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Energy 

(d^)^ 

2 

-.011214131 

- 435 

4 

. 010811795 

- 332 

Pl'fl 

. 007551306 

- 202 

(P2)' 

3 

007454889 

- 140 

(dj)^ 

5 

007170236 

- 181 

(d^)^ 

3 

006593000 

- 142 

Ssdi 

2 

005423758 

- 83 

3 

004096498 

- 72 

(d^)^ 

4 

-.003030145 

- 29 

P1P2 

2 

-. 002673776 

- 14 

^1^2 

5 

.002190747 

- 23 

(fl)^ 

5 

. 002039872 

- 26 

(d2)^ 

2 

-.001747914 

- 20 

iv^)^ 

4 

. 001610092 

- 15 

Pz^i 

1 

-.001590192 

- 10 

(fj)^ 

3 

-. 001533352 

- 16 

did2 

8 

. 001491326 

- 16 

TABLE  9 


EXPANSION 

OF  HF  INTERACTING  (s^p^) 

EXCITATIONS 

Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Energ 

=3Pz 

3 

. 054006778 

- 7130 

"3P2 

2 

. 043082303 

- 4702 

6 

. 023381234 

- 2950 

P2'^l 

6 

. 017014787 

- 912 

P3‘*l 

6 

015566722 

- 1121 

5 

015280189 

- 1314 

P2'*2 

6 

012395997 

- 957 

P2‘^l 

5 

. 010074992 

- 245 

Pjdl 

5 

-.008357891 

- 287 

P2‘^2 

5 

008274255 

- 444 

*4P3 

3 

007624900 

- 436 

*4P3 

2 

005580601 

- 212 

V2 

6 

005286835 

- 325 

"^2*2 

5 

. 004135588 

- 211 

P4'*2 

6 

. 002770989 

- 157 

S4P2 

3 

. 002725362 

114 

"3P3 

2 

. 002670829 

- 128 

Vl 

6 

001611094 

72 

= 3P4 

2 

. 001520327 

45 

TABLE  9 continued 


C onfiguration 

Degenerate 

Element 

Eigenvector 

Partial  Energy 

S3P4 

3 

001148072 

29 

P4^2 

5 

. 000947765 

24 

S4P4 

2 

000788093 

22 

S5P4 

2 

000474289 

16 

TABLE 

10 

EXPANSION 

OF  HF  NON-INTERACTING 

(s^Pj)  EXCITATIONS 

Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Energy 

4 

. 008633712 

- 297 

^2^2 

4 

. 005821690 

- 130 

difi 

3 

. 002426237 

30 

P2^2 

2 

002187318 

21 

76 


TABLE  11 


EXPANSION  OF  (p^)^  EXCITATIONS 


C onfiguration 


(P7)' 


Degenerate 

Element 


(di) 

(fl)^ 


(d^)' 

(P3)' 


Eigenvector 

-.043150565 
034767021 
-.007932459 
-. 006720458 
-. 004568322 


Partial  Energy 

- 3682 

- 3886 

- 400 

- ’ 276 

- 140 


-.001929488 
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TABLE  12 


EXPANSION  OF  L SHELL  SINGLE  EXCITATIONS 
(s^)  Excitations  . 


C onfiguration 

Degenerate 

Elencient 

Eigenvector 

Partial  En 

3 

119215608 

-28194 

"l 

2 

-.069525439 

- 6399 

"3 

1 

. 034055087 

- 1432 

1 

. 013427695 

- 148 

2 

. 009018988 

0 

'^2 

3 

-. 005997222 

536 

"2 

2 

-. 004227473 

166 

"5 

1 

-. 001661924 

11 

"a 

1 

-. 001289931 

4 

(Pj^)  Excitations. 

Configuration  Degenerate 

E lement 


Pz 


Eigenvector 


005957327 


Partial  Energy 
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TABLE  13 


EXPANSION 

OF  TRIPLE  AND 

QUADRUPLE 

L SHELL  EXCITATIONS 

2 

s^(Pj^)  Excitations  . 

Configuration 

Degenerate 
E lement 

Eigenvector 

Partial  Energy 

1 

. 006633999 

- 173 

(p/dj 

3 

. 005962582 

- 123 

2 

. 003630582 

- ■ 40 

S3P2P3 

3 

002483294 

30 

"3P2P3 

2 

-.002092872 

21 

S4(P2f 

2 

001821914 

16 

(s-jf"  Pj^  Excitations  . 

C onfiguration 

Degenerate 

Element 

Eigenvector 

Partial  Energy 

6 

004220961 

67 

(83)^  P2 

003809827 

33 

10 

. 002476778 

- 23 

S3S4P2 

2 

. 002412882 

22 

5 

002142162 

19 

P3 

. 002116743 

19 

(P2)^ 

1 

. 001848937 

13 

(P2  f"  P3 

8 

001839479 

16 

TABLE  13  continued 


(s^r  Pj^  Excitations  (continued). 

Configuration  Degenerate  Eigenvector 

Element 


1 


-.001794829 


Partial  Energy- 
17 


Excitations  . 

Configuration  Degenerate  Eigenvector  Partial  Energy 

Element 


(S3?  (P2? 


. 003620671 
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TABLE  14 


EXPANSION  OF  (s^)^  EXCITATIONS 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Er 

(P4)^ 

5 

-.010578751 

- 6996 

P3P4 

8 

-.009548171 

- 4128 

(85)^ 

-.  008621388 

- 6863 

S4S5 

2 

-.  007052000 

- -2845 

P2P4 

8 

-. 006718510 

- 1832 

P4P5 

8 

-. 005810230 

- 4225 

®3®5 

2 

. 005628016 

- 1651 

P2P3 

8 

-.005515289 

- 898 

(P3)^ 

5 

-. 005166094 

- 887 

"3®4 

2 

. 004672690 

- 683 

(d3)2 

6 

-. 004532230 

- 2554 

(P5)' 

5 

-.003895650 

- 2530 

P1P4 

3 

-.003827283 

- 573 

(34)^ 

-.003680683 

- 491 

P1P3 

3 

-.003322031 

- 304 

(P^)' 

5 

-. 003017023 

- 239 

(33)2 

-. 002982826 

- 244 

P1P2 

3 

-.002552750 

- 162 

"3^6 

2 

. 002126254 

28 

81 


TABLE  14  continued 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Ene 

P3P5 

8 

-.002058667 

- 482 

2 

001958731 

69 

(s^)2 

-.001732786 

- 411 

6 

-. 001523276 

- 650 

(^4)^ 

6 

-.001283660 

- 596 

P4P6 

8 

-.001239151 

173 

5 

-.001212527 

- 537 

2 3 

10 

-. 001072814 

97 

S5S7 

2 

. 001069407 

72 

2 

-. 001057058 

373 

P3P6 

8 

-. 000799082 

1 

P2P5 

8 

-. 000741025 

99 

(Pl)^ 

. 000722793 

16 

S4S7 

2 

. 000543591 

13 

"iS 

10 

. 000499335 

19 

(37)^ 

-.000367178 

96 

P1P5 

3 

-.000355210 

25 

P5P6 

8 

. 000299144 

- 143 

TABLE  15 


EXPANSION  OF  (s^)  EXCITATIONS 


C onfiguration 

Degenerate 

Element 

Eigenvector 

Partial  En 

3 

. 003580933 

- 183 

3 

-.002892015 

- 125 

"3 

1 

-.002339943 

72 

1 

. 001284268 

25 

d 

2 

2 

-.001216549 

22 

d 

'1 

2 

. 000577873 

20 

83 


TABLE  16 


EXPANSION  OF  (3^82)  EXCITATIONS 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Er 

P1P2 

9 

004669457 

- 312 

P1P3 

8 

-.004322516 

- 365 

P2P4 

26 

004273199 

- 533 

P1P4 

8 

-.003608779 

- 408 

(P2)^ 

7 

. 003568234 

- 179 

P2P3 

26 

-. 003448177 

- 226 

P3P4 

26 

-. 002988004 

- 263 

P2P3 

27 

. 002811528 

- 129 

(85)^ 

2 

-.002805974 

665 

P1P3 

9 

-. 002755666 

- 143 

S4S5 

6 

-. 002730652 

- 340 

s s 
3 4 

6 

. 002619509 

- 122 

(Pl)^ 

2 

-. 002603988 

92 

(P4)^ 

7 

-. 002523425 

- 353 

P1P2 

8 

-.002192258 

87 

P4P5 

27 

-. 001976463 

- 405 

(S4) 

2 

-. 001739677 

71 

®3®5 

6 

. 001660794 

- 156 

(53)^ 

2 

-. 001591106 

43 

84 


TABLE  16  continued 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Ene 

P3P4 

27 

-.001439866 

83 

(d^)^ 

7 

-.001354707 

32 

(dj)2 

7 

-. 001165515 

16 

®4®6 

6 

-.001051689 

16 

(Ps)^ 

7 

-.001034600 

- 181 

P4P5 

26 

-.000879486 

65 

P3P5 

27 

-.000795788 

58 

(d3)2 

7 

. 000689309 

53 

"5^6 

6 

-.000684705 

73 

(s^)2 

2 

-.000639710 

40 

P2P5 

27 

-. 000629025 

26 

5 

-.000547538 

15 

P4P6 

27 

-.000367892 

21 

(P6>^ 

7 

-. 000325176 

39 

(d^)^ 

7 

. 000306196 

32 

P5P6 

27 

. 000165226 

24 

85 


TABLE  17 

EXPANSION  OF  (s^p^)  EXCITATIONS 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Ene 

P d 
^3  2 

6 

. 004563477 

- 405 

P d 
^2  1 

6 

-.004451068 

- 336 

Ps^i 

6 

-.004377461 

- 395 

"3P3 

2 

. 004264900 

- 314 

P2^2 

6 

. 004106014 

- 272 

^^2 

2 

. 003652877 

- 206 

P4^^2 

6 

. 003157288 

- 321 

^P4 

2 

. 002926388 

- 251 

"4P3 

2 

-.002810999 

- 163 

"4P4 

2 

-.002582116 

- 218 

P d 

P4  1 

6 

-.002575422 

- 232 

P d 

P4  3 

6 

-. 002544174 

- 429 

"4P2 

2 

-.002257406 

89 

^SP4 

2 

-. 001978537 

- 205 

P3^3 

6 

-.001492484 

- 116 

P2^2 

5 

. 001438683 

32 

^^3 

2 

-. 001219349 

54 

®6P4 

2 

-. 001181885 

44 

P3S 

5 

-.001002019 

63 

TABLE  17  continued 


Configuration 

Degenerate 

Element 

Eigenvector 

Partial  Energy 

^2^1 

6 

. 000984567 

19 

5 

-.000932689 

52 

2 2 

6 

. 000931386 

20 

P d 
P4.  3 

5 

000899235 

58 

"6P3 

2 

-.000884990 

26 

P4^4 

6 

-.000831561 

47 

P3^^2 

5 

. 000820465 

12 

^^3^2 

6 

-. 000630841 

27 

P2^^3 

6 

-. 000612143 

25 

P3^4 

6 

-.000571633 

25 

Ps^ 

6 

-. 000552256 

55 

V3 

6 

-. 000532533 

52 

dof  0 
3 3 

5 

. 000384394 

27 

^Ps 

2 

-.000378560 

27 

^SP2 

2 

-.000305469 

8 

87 


TABLE  18 


EXPANSION  OF 

TRIPLE  AND 

QUADRUPLE  KL 

SHELL 

EXCITATIONS 

2 2 

(Sj^)  (s^)  Excitations  . 
Configuration  Degenerate 

Eigenvector 

Partial 

. Energy 

(Pl)^P4>^ 

Element 

2 

001462661 

_ 

125 

(Pl»^P3P4 

7 

001328775 

- 

76 

001039978 

- 

106 

(Pl  )^P2P4 

7 

000901142 

- 

30 

2 

000771583 

- 

40 

(Pl)?  (Pj)^ 

2 

000734083 

- 

16 

(Pi>^  P2P3 

7 

000725990 

- 

15 

(Pl^  P4P5 

7 

000688048 

- 

64 

S3S4(pj)^ 

2 

. 000661616 

- 

12 

(Pl)^d3f 

2 

000621451 

- 

45 

®3  ®5(Pl ^ 

2 

. 000595050 

- 

22 

18 

. 000553192 

- 

14 

(Pl)^P5)^ 

2 

000516001 

- 

43 

(Plf  P3P5 

7 

-.000207359 

- 

6 

(Pl)^  P5P6 

7 

. 000091817 

- 

6 

S5S^(Plf 

2 

000080621 

4 

88 


TABLE  18  continued 


Excitations  . 

C onfiguration 

Degenerate 

Element 

Eigenvector 

Partial 

Energy 

(S5)^dj 

3 

000925887 

- 

87 

(P4)^d^ 

8 

-.000824538 

- 

34 

(p/  <»i 

4 

. 000756670 

- 

28 

(P4)^  dj 

11 

000734064 

- 

27 

5 

000630004 

- 

20 

4 

. 000611896 

- 

25 

(P4f 

1 

. 000601056 

- 

18 

(85)^  d^ 

2 

000540972 

- 

30 

SjSjdi 

4 

-.000480622 

- 

14 

(si  )^S2Pi 

Excitations  . 

Configuration 

Degenerate 

Element 

Eigenvector 

Partial 

Energy 

S3P2(P4)^ 

2 

. 000624209 

- 

20 

53(85)  P2 

3 

-.000423153 

- 

18 

(Sj)?  (p,f 

Excitations  . 

C onfiguration 

Degenerate 

Element 

Eigenvector 

Partial 

Energy 

(P2>^  (P4)^ 

5 

. 000579111 

- 

17 

89 


The  Intershell  expansion  (Tables  16  and  17)  is  slowly  con- 
vergent; all  its  terms  are  HF  interacting. 

Table  18  shows  the  triple  and  quadruple  excitations  which 
involve  orbitals  of  both  shells.  The  {s^f  excitations  preserve 

the  S , S coupling  of  the  original  orbitals,  while  the  only  occurring 
(Sj^)  (Pj)  excitation  maintains  a ^ through  the  scheme  ^S(p^)? 

As  mentioned  in  Chapter  II,  it  can  be  seen  from  the  tables 
that  a close  agreement  exists  between  the  eigenvector  components  of 
the  triple  and  quadruple  excitations,  and  the  products  of  the  eigen- 
vector components  of  the  single  and  double  excitations  from  which 
they  can  be  formed. 


CHAPTER  VI 


DISCUSSION 

The  interpretation  and  utilization  of  the  results  obtained, 
involves  three  main  areas  of  interest,  namely,  (i)  the  energy, 

(ii)  the  analysis  of  the  wave  function  in  terms  of  constituents 
which  are  transferable  from  one  system  to  another,  and,  (iii)  the 
computation  of  properties  such  as  those  mentioned  in  the  Introduc- 
tion. 

We  shall  here  discuss  in  detail  only  the  first  of  these 
points.  The  second  problem,  (ii),  is  presently  being  investigated: 
in  the  next  section,  the  tranferability,  within  chemical  accuracy, 
of  K shells,  shall  be  shown.  The  systematic  investigation  of  the 
third  problem,  (iii),  concerns  not  only  the  computation  of  average 
values,  but  also  the  determination  of  how  the  latter  are  affected  by 
various  truncations  of  the  wave  function.  For  example,  since  the 
.electronic  spin  density  at  the  nucleus  depends  only  on  the  eigen- 
vector components  of  the  K shell  single  excitations  (Table  13),  one 
may  look  for  a way  to  truncate  the  total  wave  function,  which 
leaves  this  part  invariant. 


90 


91 


6.1.  The  Correlation  Energy 

It  is  very  common  practice  to  talk  of  the  energy  results  in  terms 
of  the  percentage  of  correlation  energy  (CE)  obtained.  However, 
since  the  definition  of  CE  involves  the  total  non-relativistic  energy, 
which  is  not  well-known,  its  use  as  a "yardstick"  in  which  to 
measure  the  goodness  of  wave  functions  has  to  be  made  with  some 
re  servation. 

Clementi  has  estimated  the  CE  of  the  P state  of  C to  be 

(38) 

0.158  a.u.  This  value  is  obtained  by  computing  the  average 

value  of  a simplified  relativistic  Hamiltonian,  with  respect  to  the 
HE  wave  function.  The  CE  shall  be  estimated  here  on  the  basis 
of  our  calculation  and  of  a study  of  the  errors  involved  in  the  vari- 
ous truncations. 

The  difference  between  the  CE  of  the  K shell  of  C and  that 
4't 

of  C , is  due  mainly  to  the  following  effects:  (i)  the  exclusion  of 
the  s^  orbital  from  the  expansion  of  the  K shell  of  C , which 
amounts  to  0.00257  a.u.  (see  Table  19):  (ii)  the  partial  exclusion 
of  the  pj  orbital  from  the  K shell  expansion  of  C , due  to  the  pres- 
ence of  the  (Pj^)^  pair  in  the  L shell  (about  0.00059  a.u.);  and  (iii) 
the  absence,  from  the  C^  K shell  expansion,  of  the  s^  single 
excitations  which  appear  in  the  corresponding  C wave  function,  due 
to  the  fact  that  they  do  not  preserve  the  coupling  of  the  K shell. 
This  energy  effect  (about  0.00045  a.u.)  opposes  those  of  (i)  and  (ii). 


ANGULAR  ENERGY  LIMITS  FOR 
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Assuming,  then,  that  all  the  difference  between  the  CE's 

4+ 

of  the  K shells  of  C and  C can  be  accounted  for,  their  errors 
should  be  the  same:  the  total  error  in  is  obtained  in  Table  19 

on  the  basis  of  the  angular  energy  limits. 

The  calculation  of  the  error  in  the  L shell  CE  is  done  in  an 
entirely  different  manner.  The  error  due  to  the  truncation  of  the 
L shell  Cl  expansion  is  evaluated  by  adding  the  partial  energies  of 
the  terms  which  were  dropped  during  the  search.  The  effect  of  g-, 
h-,  and  higher  symmetry  orbitals  is  estimated  for  each  electron 
pair,  according  to  an  empirical  procedure  developed  by  Bunge, 
which  gives  very  accurate  predictions  for  the  case  of  the  He  iso- 
electronic  series.  It  is  illustrated  in  Table  20  for  the  (s  exci- 
tations.  It  is  based  on  the  existence  of  a general  pattern  for  the 
contributions  e(i,£)  of  the  different  orbitals  (i,  f)  to  the  energy,  when 
these  are  written  in  the  manner  shown  in  Table  20.  The  energy 
e(i,£)  is  obtained  by  adding  the  partial  energies  of  the  excitations 
in  which  (i,  £)  occurs;  only  those  excitations  of  (i,  f)  are  considered, 
in  which  (i,  f)  occurs  paired  with  itself  or  with  another  orbital  (i',£') 
for  which  e(i',f')>  e(i,£)  . One  finds  that: 

a)  e(i,  £ + 1)  < e(i,  f) 

b)  e(i,£  + l)  >e(i  + l,f) 

c)  the  vertical  ratios  e(i,  f)/e(i  + l,f)  usually  decrease  smooth- 
ly for  a given  i value,  for  increasing  £ . 

Once  a few  of  the  e(i,  f)'s  have  been  obtained,  the  remaining 
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TABLE  20 

DIAGRAM  TO  ESTIMATE  e(i,  !)  VALUES  FOR  (s^)  EXCITATIONS 

d =4100^ 

1 

d^=  580  f =1000 
2 1 

ds=  20  1^=  150  g^=  300 

d^=  0 f^=  10  g^=  70  h^=  180 


a.  In  units  of  10  a.  u.  (C). 
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ones  are  predicted  according  to  the  rules  indicated  above:  in  this 

way,  we  have  estimated  the  effect  of  the  g-,  h-,  and  higher  orbit- 

als  on  the  expansion  of  (s^)'^,  to  be  of  the  order  of  0.0010  a.u. 

Analogously,  we  find  0.00045  a.u.  and  0.0030  a.u.  for  the  pairs 

2 

(Pj^)  and  S2Pj^  , respectively. 

For  the  Inter  shell,  similar  arguments  lead  to  an  estimate 
of  the  error  in  the  CE,  of  0.00020  a.u. 

Finally,  all  the  errors  are  combined  as  described  in  Table 
21 , to  give  CE=-0.1546±  0.0010  a.u. 

6.2.  Comparison  with  other  Works 

As  compared  with  the  enormotis  amount  of  literature  on  the 

ground  state  of  Be,  there  is  indeed  very  little  published  on  the 

ground  state  of  C.  In  Table  22,  we  have  given  the  results  for 
3 

C(  P),of  variational  calculations  which  go  beyond  the  HF  approx- 
imation. Since  some  of  them  are  concerned  only  with  the  L shell 
Cl  expansion,  the  energy  results  for  both  the  latter  and  the  total 
Cl  wave  function  are  displayed. 

Boys's  calculation  is  the  first  extensive  Cl  w'ork  on  C. 

Weiss'  emphasis  is  related  to  the  computation  of  transition  proba- 
bilities,  ' ' rather  than  the  determination  of  wave  functions:  most 

of  the  im±provernent  of  our  L shell  expansion  over  his,  comes  from 
the  partition  of  degenerate  spaces.  Bagus  and  Moser  perform, 

a multiconfiguration  self-consistent  field  calculation,  aiming  at  the 
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TABLE  21 

ESTIMATE  OF  THE  TOTAL  NON -RELATIVISTIC  ENERGY 


Ene  rgy 

Full  Cl  estimate^ 

-37. 83620 

4+ 

K shell  error  = K shell  error  C 

- .00165 

L shell  and  Inter  shell  spdf  saturation^ 
remainders 

- .00067 

Contributions  of  g-,  h-,  and  higher 
symmetry  orbitals  to  (82)^  excitations 

- .00100 

2 

Same,  to  (p^ ) excitations 

- . 00045 

Same,  to  (s2Pj^)  excitations 

- .00300 

Same,  to  Intershell  excitations 

- .00020 

Total  non- relativistic  energy  (estimate) 

-37. 84317 

HF  energy*^ 

-37. 68861 

Correlation  energy  (estimate) 

- .1546± 

a.  See  Table  6. 

b.  See  Table  19. 

c.  See  Section  5.  2. 

d.  Reference  4. 


97 


(M 

H 

< 

H 


W 

H 

< 

Eh 

W 

Q 

is 

o 

o 

u 

h 

o 

w 

o 

l-H 

H 

< 

u 

C 

u 


u 


u 


o 


o 

C 

!3 


(U 

> 

nJ 


4-> 

0 

H 


>N 

W) 

h 

0) 

ti 

H 


01 

a 

h 

(U 

Eh 


CO 

•H 

CO 

PQ 


't 

r^ 

ro 

I 


m 


XI 


a 

CO 

DO 


00 

r- 

co 

CO 

00 

CO 

I 


CO 

(M 


CO 

T3 

(i 

vD 

CO 


>s 

&0 

<1) 

!h 

H 


ID 

00 

o 

00 

00 

fO 

vO 

00 

(N3 

in 

r^ 

00 

00 

r- 

CO 

CO 

CO 

CO 

I 1.1 


U 


lU 

W 


(U 

H 


CO 

cci 

m 


o o o 

CO 


»4H 

HH 

(VJ 

(NJ 

(NJ 

•k 

•k 

•k 

X) 

'V 

(VJ 

ro 

(NJ 

•k 

•k 

•k 

Oh 

fX 

a 

Ph 

(VJ 

CO 

•. 

•k 

•k 

•k 

CO 

CO 

CO 

CO 

(NJ 

lO 

c^ 

0) 

CO 

O 

t3 

!P 

ccj 


C 

O ^ 

•X3  C 

.2 


P! 

U 


<0 

(0 

pS 

CO 

CO 

CO  kii 

bO 

>> 

•rt 

(t 

0 

CQ 

CQ 

Eh 

C 

o 


P3  -<-> 

U CO 

Xi  l> 

cfl  ^ 


Eh 


CD 

CT' 

. 

CO 

00 

lU 

0 

4) 

u 

CJ 

CJ 

c 

p: 

C 

0) 

V 

CD 

u 

(U 

<u 

0) 

v< 

m 

(D 

cu 

<1> 

Pi 

Pi 

Pi 

c3 

u 

98 

obtaining  of  a compact  representation  of  a Cl  wave  function  which 
is  restricted  both  in  the  type  of  configurations  involved  and  in  the 
internal  couplings  allowed. 

In  Table  23,  we  have  given  the  partition  of  the  calculated  CE 
in  terms  of  the  e parameters  of  HF  orbitals  or  of  pairs  of  HF  orbit- 
als: these  are  defined  as  the  sum  of  the  partial  energies  of  the 
corresponding  excitations,  for  example  e(s2.p^)=S  where  the 

summation  goes  over  all  the  excitations. 

McKoy  and  Sinanoglu  have  made  estimates  of  pair  corre- 
lation energies,  using  the  many-electron  theory.  Their  results 

correspond  to  a different  partition  of  the  CE,  but  some  fragments 
can,  however,  be  compared  with  ours:  in  particular,  their  e(2s^) 

= -.0150  a.u.,  is  equivalent  to  our  contribution  from  the  configuration 
2 4 

(Pj ) • Also,  their  combined  pair  correlation  energies  s{2p— * 

2,2  2 

Is  2s  )+  e(2p  ) , which  amount  to  -.0851  a.u.,  is  closely  related 
to  our  e(s2,P]^)+  e(s  ^ , pj^ ) + e(s2)  + e (p^  , p^ ) which  gives  -.0655  a.u. 

6.3.  Conclusions 

3 

The  results  obtained  for  the  P state  of  the  carbon  atom 
have  surpassed  even  our  most  optimistic  expectations.  Although 
different  states  of  different  atoms  may  present  new  difficulties, 
several  problems  have  been  approached  here,  which  are  of  a general 
nature. 

The  three  main  aspects  of  our  calculation  are:  (i)  the  opti  - 
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TABLE  23 

PARTITION  OF  THE  CALCULATED  CORRELATION  ENERGY 


a 

Energy 

^HF 

-37. 688547 

£(s2» 

- .028626 

e(s2.  P^) 

- .021625 

e(Pl.Pl) 

- .008433 

e(s^) 

- .035456 

s(sj  , ) 

- .040575 

e(3j) 

- . 000447 

e(s^,  s^) 

- .005446 

e(sj , p^) 

- . 004623 

E(total) 

-37. 833778 

a.  Reference  25. 
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mization  of  the  STO  basis,  together  with  the  evaluation  of  the  satu- 
ration remainders;  (ii)  the  simplification  of  the  description  of  the 
full  Cl  expansion  by  means  of  NRO's;  (iii)  the  optimum  partition  of 
the  degenerate  spaces,  which  is  especially  significant  for  open  shell 
systems. 

Because  the  relativistic  corrections  have  not  been  considered 

in  this  work,  a comparison  between  our  energy  results  and  the 

"experimental”  energy  E is  of  doubtful  significance.  The  latter 

can  be  found  by  adding  the  "experimental"  values  for  the  successive 

(42) 

ionization  potentials  of  carbon,  except  for  the  ionization  potential 

4+ 

of  C , which  is  given  more  accurately  by  the  theoretical  work  of 
(43) 

Pekeris  : we  find  E =- 37.  85638±  0.  00002  a.u.  (C). 

exp. 

With  a large  electronic  computer,  it  seems  possible  to  obtain 
results  of  comparable  accuracy,  in  short  periods  of  time.  The 
computation  of  the  1-matrix  of  our  wave  function  is  under  way.  We 
plan  to  calculate  next  the  other  two  states  of  C arising  from  the 
HE  configuration  (s^)^(s2)^(p^)^,  namely,  the  lowest  ^ and  states; 
the  corresponding  wave  functions  will  be  used  to  obtain  relativistic 
corrections  to  the  energy,  hyperfine  constants,  and  other  properties. 
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